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SLENDER-BODY THEORY BASED ON APPROXIMATE SOLUTION 
OF THE TRANSONIC FLOW EQUATION 

By John R. Sprf.iter and Alberta Y. Alksne 


SUMMARY 

Approximate solutions of the nonlinear equations 
of the small disturbance theory of transonic flow are 
found for the pressure distribution on pointed slender 
bodies of revolution for flows with free-str earn Mach 
number 1 , and for flows that are either purely sub- 
sonic or purely supersonic . These results are ob- 
tained by application of a method based on local 
linearization that was introduced recently in the 
analysis of similar problems in two-dimensional 
flows. The theory is developed for bodies of arbitrary 
shape , and specific results are given for cone-cylinders 
and for parabolic-arc bodies at zero angle of attack . 
All results are compared either with existing theo- 
retical results or with experimental data . 

INTRODUCTION 

This paper is concerned with the prediction of 
the pressure distribution on slender pointed bodies 
of revolution in flight at zero angle of attack at 
high subsonic Mach numbers, at low supersonic 
Mach numbers, and at Mach numbers near 1. 
The solution of such problems is of interest not 
only because of the frequent use of bodies of revolu- 
tion in practical applications but also because 
knowledge of the aerodynamic properties of a 
body of revolution in axisymmetric flow taken to- 
gether with the transonic area rule and equivalence 
rule described in references 1, 2, 3, and elsewhere 
permits the ready calculation of the aerodynamic 
properties of a wide class of wings, bodies, wing- 
body combinations, etc., having the same longi- 
tudinal distribution of cross-section area as the 
body of revolution. 

The analysis is based throughout on the small 
disturbance theory of transonic flow. This the- 
ory, at least as formulated in this paper and in 


several earlier papers of the present authors, is not 
restrictive in its range of applicability as might he 
inferred from its name, but is equally applicable 
to subsonic and supersonic flows. It is, moreover, 
the simplest theory proposed to date that is 
capable of yielding reliable results throughout this 
entire speed range. Difficulties arise in the solu-. 
tion of the equations of transonic flow theory, 
however, because the governing partial differential 
equation is nonlinear and of mixed elliptic-hyper- 
bolic type if regions of subsonic and supersonic 
flow both appear in a single flow field. Inasmuch 
as the mathematical theory of such equations is 
still in a rather early state of development and 
methods have not yet been discovered for the 
exact solution of such equations, the engineer is 
forced to turn either to experiment or to approxi- 
mate methods for the solution of practical prob- 
lems. Of the various methods that have been 
proposed for the approximate solution of the equa- 
tions of transonic flow theory, the most successful 
and also the most versatile method is that de- 
scribed recently for two-dimensional flows in refer- 
ence 4. This was demonstrated in reference 4 by 
the presentation of results for a large number of 
specific applications for which experimental and 
other theoretical results were available for com- 
parison. Additional results for both two-dimen- 
sional and axisymmetric flows with free-stream 
Mach number 1 are given in reference 5. The 
latter account is brief, however, and it is the pur- 
pose of the present paper to provide a more de- 
tailed account of the extension and application to 
axisymmetric flows of the method described origi- 
nally in reference 4. Inasmuch as this method of 
approximation is of somewhat novel character and 
the mathematical basis for the procedure is only 
incompletely understood, the results are again 

l 
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evaluated by extensive comparisons with experi- 
mental results and also with other theoretical re- 
sults whenever available. 
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PRINCIPAL SYMBOLS 


parameter defined in equation (50) 
Euler’s constant ~ 0.5772 15665 . . . 

pressure drag coefficient, ( p j%) J! 2 S 


pressure coefficient, 


P~P - 
(pJ2) U 


2 


pressure drag 

maximum diameter of body of revolution 
interval between successive values of x in 


numerical integration 

MJ (t+1) 


length of body 
local Mach number 
free-stream Mach number 
local static pressure 
free-stream static pressure 
radius of body 

area of body cross section, ttR 2 
free-stream velocity 

perturbation velocity components parallel 
to x and r axes, respectively 
value of u obtained by solution of equation 
( 16 ) 

value of u obtained by solution of equation 
(35) 

value of u obtained by solution of equation 


(47) 

cylindrical coordinates where x extends in 
tiie direction of the free-stream velocity 
ratio of specific heats, 1 .4 for air 
difference between ‘‘predictor” and “cor- 
rector” for numerical integration (see 
eqs. (23) and (24)) 
l-MJ-ku 
MJ-~\+ku 



free-stream density of air 
thickness ratio, j 
perturbation velocity potent ial 


SUBSCRIPTS 

values associated with critical Mach num- 
ber 


i values associated with incompressible flow 
or with 0 

m values associated with maximum diameter 
of body 

x,r differentiation with respect to the variable 
in question 

FUNDAMENTAL EQUATIONS AND BOUNDARY 
CONDITIONS 

STATEMENT OF PROBLEM 

Consider the steady axisymmetric flow of an 
inviscid compressible gas past an arbitrary slender 
pointed body of revolution, and introduce cylindri- 
cal coordinates x and r with the x axis parallel to 
the direction of the free stream, as illustrated in 


A /-/?(*) 



Figure 1. -View of body and coordinate system. 


figure 1. Let the free-stream velocity and density 
be f/oo and p*,, the perturbation potential be *p, 
and the perturbation velocity components parallel 
to the x and r axes be <p x or u , and <p r or p, where 
the subscript indicates differentiation . The bound- 
ary conditions require that the perturbation veloci- 
ties vanish at infinity, and that the flow be tangential 
to the body surface. The first condition indicates 
that ip is constant at infinity. The second condi- 
tion can be approximated for smooth slender 


bodies by 


(r<pr) T -Q=U„Ii 


dRJJ a 
dx 2 tt 


dS 

dx 


( 1 ) 


where R and £ represent, respectively, the ordi- 
nates and cross-section area of the body. The 
relation between the pressure coefficient C p and 
the perturbation velocity components is likewise 
approximated and reduces to the following form 
for points near the body surface: 


C, 



o <Px <Pr 2 

U m UJ 


( 2 ) 


The corresponding expression for points on the 
bod) 7 surface is 



cr 


(3) 
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These relations are familiar from linear theory but 
apply equally for transonic theory. The differ- 
ential equation for <p is not the same as in linear 
theory, however, but is 


/I T f 2\ | 1 I 7 "hi) 7. 

)<p xz - (Pr-r<prr= <Px<Pzx = k<Pz<Pxz 


( 4 ) 


where M m is the Mach number of the undisturbed 
flow and y is the ratio of specific heats (1.4 for air) . 
It is useful to note that the total coefficient of p XT} 
that is, 1 —M a 2 —k<p X) corresponds in the present 
approximation to 1 —M 2 where M is the local 
Mach number. 

In transonic and supersonic flows, it is also 
necessary, in general, to provide appropriate 
relations for the discontinuous changes in velocity 
that occur at shock surfaces. The necessary 
equations, when simplified to the form consistent 
with the approximations of transonic flow theory, 
reduce to 


(1 —MJ) (u a -u b y+ (v a ~~v b ) 2 ' 

= /w a +wA 

L ~ \ z y (5) 

<Pa~<Pb 
u a >u b 

where the subscripts a and b refer to the values 
immediately upstream and downstream of the 
shock surfaces. 

Solutions for problems of transonic flow theory 
must satisfy not only equations (1) through (5), 
but certain additional requirements concerning 
the regions of influence and dependence associated 
with subsonic and supersonic flows. Although a 
detailed discussion of these conditions is com- 
plicated by the presence of both subsonic and 
supersonic flows in a single flow field, the principal 
conclusion is that the usual conditions apply 
locally and must be satisfied in order to obtain 
solutions that are physically realistic (see ref. 6, 
pp. 448-453, for an elementary discussion of this 
mat ter) . 

SIMILARITY RULES 

The equations for axisymmetric transonic flow 
enumerated above contain a similarity rule that 
associates the pressure fields of affinely related 
families of bodies. Various members of such a 
family of bodies may be of different length l and 


thickness ratio t, but must have ordinates given 
by an expression of the form 

*(?)-<?) < 6) 

where /is the same function of xjl for all members 
of a single family. If the subscripts 1 and 2 refer 
to different members of a given family, the simi- 
larity rule states that the pressure coefficients C P 
at corresponding points defined by 

J "2 r-i /M», 2 ( Ti + 1) r,r, 

1 2 W h 1W, W 

are related according to 

( 8 ) 

provided the free-stream Mach numbers of the 
two flows are in accordance with the following 
expression : 


i Af.,* /T 2 y i | 2 
A/» a 2 (7 2 +l) W A/^Cy. + i) ■ ' 

The similarity rule given above cannot be used 
directly to relate the surface pressures on bodies 
having different thickness ratios because the 
ordinates of related bodies given by equation (6) 
do not conform with the relationship for cor- 
responding points given by equation (7). Thus 
the r coordinate of a point in the vicinity of body 
2 that corresponds to a point on the surface of 
body 1 is given by 


r 2 _ y«i 2 (7i+i) / r,v 

n 2 va/„ 3 2 (7 2 +i) w 


( 10 ) 


Oswatitsch and Bemdt have shown (ref. 7) 
that a similarity rule can be established for the 
surface pressures on affinely related bodies of 
revolution if it is assumed that the longitudinal 
perturbation velocity component u is given by 
an expression of the form 

'+»<*> < l » 

in the vicinity of the body. This relation permits 
the calculation of the difference in pressure 
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between the point r 2 and the surface of body 2. 
The similarity rule for the surface pressures follows 
immediately. It. states that 



provided still that the free-stream Mach numbers 
are in accordance with equation (9). A similarity 
rule for drag D can be obtained by combination 
of equation (12) and the following expression for 
the drag coefficient C D : 


C D = 


T) 


P co 


:-s„. 


r 

n J() 


n <!S , 
C„ -t- dr 


~S m J o ' ” dr ' 


(13) 


It is 


Cr 



, sni) | n 

1 2irS mi TjS'JL a *(V2+l) 


(14) 


where the prime indicates differentiation with 
respect to x and S m refers to the maximum cross- 
section area. This relation simplifies if the body 
is either pointed or cylindrical at the base, since 
then S'(l) is zero. 

It should he observed that the relation given 
by equation (11) follows as a direct consequence 
of the assumption implicit in the use of the 
slender-body boundary condition given by equa- 
tion (1) that at any given x the product r4> T is 
constant between the x axis and points in the 
vicinity of the body surface. The error intro- 
duced into the analysis by the use of equation 
(11) is thus of the same order of magnitude as 
that already present in the fundamental relations 
of the small disturbance theory of transonic flow 
enumerated in equations (1) through (5). It is 
important to recognize, however, that equations 
(11), (12), and (14) must be regarded as approxi- 
mate relations within the framework of analysis 
specified by equations (1) through (5), and that 
exact solutions of the latter equations will not 
obey the similarity rules perfectly. This situa- 
tion may be contrasted with that for the similarity 
rule for airfoils or for wings of finite span for 
which the similarity rules are exact relations 


within the framework of the small disturbance 
theory for transonic flow. 

A much stronger limitation that affects the 
general applicability of the similarity rules for 
axisymrnctric flows, and to which attention does 
not appear to have been directed previously, is 
concerned with the regions of influence and 
dependence. This limitation arises because the 
extent and even the existence of a region of 
supersonic flow imbedded in a subsonic flow, or 
of a subsonic region imbedded in a supersonic 
flow, are not determined by the value of the simi- 
larity parameter {MJ— 1)/[A/« 2 (y+ 1)t 2 ] alone, 
as in flows around thin wings, but by the values 
for the Mach number and thickness ratio indi- 
vidually. These effects are small at Mach 
numbers near unity for which the sonic line is 
inclined at a steep angle to the body surface, 
but become large as the inclination of the. sonic 
line approaches that of the body surface at Mach 
numbers near the lower and upper critical. The 
latter Mach numbers are defined, respectively, 
as the lowest subsonic Mach number and the 
highest supersonic Mach number at which sonic 
vcloci ty (he., uj Z r Qa — (\—M M 2 ) + 1 )]) oc- 
curs at the surface of the body. Consider, as 
an example of a case for which these effects are 
large, that results are available for a certain body 
at the lower critical Mach number and that it 
is desired to use the similarity rule to determine 
the flow around an affinely related body that is 
more slender than the first body. The flow past 
the second body will be found to contain a super- 
sonic region, but since there are no shock waves 
in the flow around the first body there will be no 
shock waves indicated in the flow around the sec- 
ond body. If, on the other hand, results are 
available for a body at a Mach number somewhat 
greater than the lower critical, and the flow con- 
tains a shock wave, application of the similarity 
rule to compute the flow around an affinely 
related body of sufficiently greater thickness that 
the flow is purely subsonic leads to the completely 
unacceptable result that there are discontinuities 
in a purely subsonic flow. It is evident that 
equally unacceptable results are indicated b} r the 
similarity rule for drag, and that related diffi- 
culties with both the details of the flow field 
and the drag are encountered at Mach numbers 
in the vicinity of the upper critical. 

The properties of the similarity rule discussed 
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in the preceding paragraph require the imposi- 
tion, for free-stream Mach numbers other than 1, 
of certain rather definite restrictions on the values 
for r, or for for which the surface pressures 
and drag are related by equations (12) and (14). 
It is clear such a restriction is necessary because 
the similarity rules contain no mechanism for the 
subtraction or addition of shock waves and the 
associated contributions to drag, as there should 
be in any proper relationship between subcritical 
and supercritical flows. These and related points 
will be examined further in the course of the 
following discussion. 

The remainder of this paper is concerned with 
the presentation of a method for the approximate 
solution of the problem described by equations (1) 
through (5), application to a number of specific 
bodies, and comparison with existing theoretical 
and experimental results. Purely subsonic flows 
are discussed first, purely supersonic flows next, 
and flows with free-stream Mach number near 
unity are treated last. 


follows immediately by differentiation. It can be 
approximated for points on the surface of a smooth 
slender body by 


TT a 


S"(r) Tn \bS 1 f 1 S' (7) -S' (i) , 

4t Atx{1—x)' Atj 0 jr — £[ 


S"(x) 

Air 


In 

r oo 


(18) 


where the subscript i refers to the values for in- 
compressible flow, or for Differentiation 

yields 


ax 4tt ax 


If, now, 1— A/a , 2 — ku is restored in place of \ E so 
that, in effect, the local value for \ E is used at each 
point, and the subscript E on u is dropped, equa- 
tion (19) becomes 


d(ulUJ _ S'"(x) 
dx Air 


In (1 —M^—ku) + 


w 

dx 


(20) 


SUBSONIC FLOWS 

APPROXIMATE SOLUTION OF EQUATIONS 

It is convenient in the analysis of subsonic 
flows to introduce the sjunbol \ E as an abbrevi- 
ation for the coefficient of <p IX 

\ E = 1 -MJ-MJ y £-<p x = 1 0 ( 15 ) 

C CO 

and rewrite equation (4) in the form 

+ (16) 


It is now assumed that \ E is neither zero nor 
infinite and that it varies sufficiently slowly that 
its derivatives can be disregarded so that it can be 
considered, temporarily, as a constant. At this 
stage, the problem is equivalent to that encoun- 
tered in linearized theory of subsonic flow around 
slender bodies of revolution (it is identical if is 
replaced by 1-M W 2 ). The solution for <p is thus 
well known (see, e.g., ref. 8, p. 188) and is given by 


u.r ^ mi 

4?r Jo \ (r— f) 2 -j-X^r 2 


(17) 


where the subscript E denotes that the solution 
refers to equation (16) of elliptic type with X^ 
constant. The corresponding expression for u E 


Equation (20) is a nonlinear ordinary differential 
equation for a on the body surface. This equation 
is of first order and the necessary value for the 
constant of int egration is determined by assuming 
that u equals u E at the point on the forebody 
where S"(x) vanishes, since u B is indicated by 
equation (18) to be independent of \ E and hence 
of Mach number at this point. Thus 

u = Ui where S"(x) = Q (21) 

The calculation of u on the surface of a body of 
revolution of specified shape is thus reduced to 
the solution of an equation of the form 

S-FM ( 22 ) 

which can be solved, at least approximately, by 
application of any of a number of standard 
techniques. The corresponding values for the 
surface pressures follow directly upon insertion of 
the results into equation (3). 

The particular method used to calculate the 
results presented in the following section of this 
report is that of Milne (ref. 9, p. 135). It is a 
step-by-step process in which the values of the 
dependent variable u are calculated one after the 
other for a sequence of equally spaced values of 
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the independent variable x. The successive 
values of x are denoted by x 0 , x u x 2f . . the 
interval is denoted by h, and the corresponding 
values of u are denoted by ?/ 0 , u l} u 2) . . . - Dif- 
ferentiation with respect to x is indicated by a 
prime. The actual calculations proceed by alter- 
nate use of the formula 

u n + 1 — u n -3 -\-^h(2u n ' — u n _ 1 * + 2 u n _ 2 1 ) ( 23 ) 

as a “predictor” and Simpson's rule 

■ h»+i — ( 24 ) 


where r represents the thickness ration djL The 
variation of cross-section area S with x is thus 

s=tR ^ AwtH 2 |j_Qy] ( 26 ) 

and the points at which S"(x) vanishes are located 

(27) 

The expression for w at zero Mach number can 
be determined directly upon substitution of the 
above relation for S(x) into equation (18) and 
carrying out the indicated operations. The result 
can be expressed simply as 


as a “corrector.” Both expressions are approxi- 
mate, but the error in the result given by the 
corrector is roughly equal to A/29 when A is the 
difference between the results given by the two 
formulas. It is assumed that the value given by 
equation (24) is correct so long as A/29 is not 
significant. Abrupt fluctuations in the values 
of A indicate the presence of errors and the calcula- 
tions should be checked. If the error A/29 
proves to be larger than the desired accuracy 
permits, it is necessary to shorten the interval h. 
Cutting the interval in half will divide the error 
by about 32. 

The computation is started at the point where 
S"(x) is zero and requires four consecutive known 
values for u. One of these is given by equation 
(21), the remaining three are determined by use 
of the first six terms of a Taylor's series for u about 
the point where S"(x) vanishes. This requires 
evaluation of the first five derivatives of u with 
respect to x at the point where S"(x) vanishes. 
The first of these is calculated directly by use of 
equations (20) and (21). The remaining four 
are calculated similarly by use of corresponding 
expressions derived from equation (20) by succes- 
sive differentiation with respect to x. 

APPLICATION TO PARABOLIC-ARC BODY 

As an application of the foregoing results, 
consider subsonic flow past a parabolic-arc body 
of revolution of length l and maximum diameter 
d. The ordinates 7? of this body are given by 

fl = 2w [?-(0’] <25 > 


where 7?(x) and S(x) now refer to the radius and 
cross-section area of a parabolic-arc body of 
revolution, and are defined as functions of x by 
equations (25) and (26). Substitution of equation 
(28) into equation (20) leads to the following 
ordinary nonlinear differential equation of first 
order for u on the surface of a parabolic-arc body 
of revolution in subsonic flow. 


d (u/UJ) _S'" 
dx 4ir 


In (1 -MJ-ku) 





The constant of integration is evaluated by im- 
posing the condition derived from equations (21), 
(27), and (28) that 


u-=0 



(30) 


The problem described by equations (26), (29), 
and (30) has been solved numerically for several 
suberitieal Mach numbers for a body of fineness 
ratio Ijd of 10. The resulting values of u/U m are 
tabulated in table I together with the values for 
the error term A/29. Results are given for only 
the front half of the body since the solution is 
symmetric around the center of the body, and the 
calculated values at x and l—x should be identical. 
It can be seen from inspection of the values for 
A/29 that the size of the interval h, which was 
maintained at 2 percent of the body length 
throughout all of the calculations, is sufficiently 
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small that inaccuracies incurred in the numerical 
solution are negligible. 

The corresponding values for C v on the body 
surface can be calculated directly from the values 
for u/U*, by means of equation (3). The resulting 
values are tabulated in table I, and are shown 
graphically in figure 2. The corresponding re- 
sults indicated by linearized compressible flow 
theory arc shown in figure 3. 

The latter results are calculated by means of 
the following expression derived from equations 
(3), (18), and (25) with \ E replaced by 1 —MJ. 


C =-^Y; 

2 A 


3+ln 


VI -MJrB 

21 


j- R ' 2 


(31) 



-L_- J J ' — 1 

.1 .2 .3 .4 .5 

x/l 


Figure 2. — Pressure distributions on a parabolic-arc body 
of revolution of fineness ratio 10 at subcritical Mach 
numbers as indicated by present theory. 

2 


It can be seen that the values for C P obtained 
by application of the present theory are qualita- 
tively similar to those indicated by linear theory, 
although somewhat more negative at all points. 
This trend is similar to that shown in reference 4 
to result upon application of the same procedures 
to two-dimensional subsonic flows. 

The results presented in table I and figures 2 
and 3 are repeated in figure 4 together with the 
corresponding results indicated by the second- 
order approximation to the solution of the equa- 
tions of transonic flow theory for flow past a 
slender parabolic-arc body of revolution. The 



Figure 3. — Pressure distributions on a parabolic-arc body 
of revolution of fineness ratio 10 at subcritical Mach 
numbers as indicated by linear theory. 


501 377 — 59 - 



8 


TECHNICAL REPORT R-2— NATIONAL AERONAUTICS AND SPACE ADMINISTRATION 



(a) A/*, =0.70 (b) M„=0.S0 (c) 0.85 (d) Af„=0.90 

(e) A/ O) =0.92 (f) M a =0.93 (g) M a =0.934 (li) 0.935 

Figure 4. — Pressure distributions on a parabolic-arc body of revolution of fineness ratio 10 at subcritical Mach numbers 
as indicated by present theory, by other theoretical results, and bv experiment. 


latter results are calculated by use of the expres- 
sion 


C =C 

V -'Z>2 P 1 


MJh I 0 

1 -MJ 




(32) 


which is obtained from a corresponding result 
given by Van Dyke in reference 10. The ex- 
pression given in reference 10 differs from equation 


(32) in that the quantity MJ { 7 + 1) is replaced 
by 7 +I, This change, which is associated with 
a corresponding difference in the coefficient k 
of the nonlinear term of equation (4), has very 
little effect on the values for C p calculated by 
means of equation (32) since the quantity involved 
is absent in the first-order term and first appears 
in the much smaller contribution of the second- 
order term. The consequences of sucli a change 
in the quantity represented by k become of con- 
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siderably greater significance in the calculation of 
transonic flows and in the determination of the 
critical Mach number, however, because the mag- 
nitude of u/U a , associated with the occurrence of 
sonic velocity is inversely proportional to Jc. As 
an example, the critical value for the pressure 
coefficient C v associated with the frcc-stream 
Mach number of 0.85 and zero <p r is —0.320 if 
MJ(y + l) is used for kU m , and —0.231 if 7+1 is 
used. Since the exact value indicated by iscn- 
tropic flow theory is —0.302, it follows that the 
formulation of transonic flow theory in which 
kU m is equal to Af„ 2 ( 7 + 1 ) is about four times as 
accurate in this respect as the formulation in 
which kU„ is equal to 7 +I. 

Experimental data from reference 1 1 by Taylor 
and McDevitt are also included in this figure for 
all subcritical Mach numbers for which results 
were measured. The tests were conducted in the 
Ames 14-foot transonic wind tunnel with a model 
that consists of approximately the forward 85 per- 
cent of a parabolic-arc body of revolution of fine- 
ness ratio 10 and length 6.67 feet. The model 
was supported from the rear by a circular cylinder 
having a diameter somewhat smaller than that of 
the base of the body. No corrections have been 
applied to the data to take account of either wind- 
tunnel wall interference or the abbreviated length 
of the body. 

It can be seen from the foregoing comparisons 
that the principal differences between the pressure 
distribution indicated by linear theory and by 
higher approximations occur near the middle of 
the body at Mach numbers near the critical. In- 
sight into the magnitude of these effects can be 
had by examination of figure 5 which shows the 
variation of C p with M „ at the midpoint of a 



Finmn 5. — Variation of pressure coefficient with Mach 
number at the midpoint of a parabolic-arc body as in- 
dicated by present theory, by experiment, and by two 
different methods of successive approximation to the 
transonic equation. 

parabolic-arc body of revolution of fineness ratio 
1 0 . Also included in this figure is a lino indicating 
the variation of the critical pressure coefficient 
associated with the occurrence of sonic velocity at 
a point, such as the midpoint of a parabolic-are 
body, at which dTtjdx is zero. The present results 
are included in this figure together with experi- 
mental results from reference 1 1 , theoretical re- 
sults indicated by equations (31) and (32), and 
by the first two steps of an alternative method of 
successive approximation described in the appen- 
dix of reference 12 that involves the solution of 
quadratic rather than linear equations at each step 
of the iteration process. The latter values for C v 
on the surface of a body of revolution are calcu- 
lated by application of the following expression 
in which the subscript N refers to the results indi- 
cated by the quadratic method and n refers to the 
nth approximation using the classical method of 
successive approximation : 




/(IR \ 2 
\dx ) 


- 2(1 -MJ) 
MJ( 7 + 1 )' 


Hi 


14 




\c v 


.+(i)’]} 


MJ(y+ 1 ) 

1 — Af 2 


(P*n 


(33) 


in which 



Values for C Pn=1 and C Pn=2 can he calculated for 
parabolic-arc bodies by application of equations 
(31) and (32). It is evident from the form of 
equation (33) that the values for C Vn approach 
those for G v in the limit of an infinite number of 
iteration steps, provided, of course, that the series 


expression for C Pn converges. If only a finite num- 
ber of terms are considered, as must be done in 
any practical application of the method of succes- 
sive approximation, the approximate values for 
the pressure indicated by C Vn and C Pw are different. 
The utility of considering the two sets of results 
simultaneously resides in the fact that the values 
indicated by C Pn tend to underestimate the true 
variation of the peak negative pressure with Mach 
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number, whereas those indicated by C Py tend to 
overestimate the true variation. The true varia- 
tion is thus bounded within well defined limits. 
It can be seen from examination of figure 5 that 
the results of the present numerical calculations 
fall within or very close to these limits. Similar 
applications of t-lie quadratic method of successive 
approximations to two-dimensional flows can be 
found in references 4 and 13. 

The problem described by equations (26) , (29), 
and (30) has also been solved numerically for sev- 
eral subcritical Mach numbers for a parabolic-arc 
body of revolution of fineness ratio 6. The result- 
ing values for u/U A/29, and C p are tabulated in 
Table II. The pressure distributions are also illus- 
trated graphically in figure 6. The corresponding 
results indicated by linearized compressible flow 



Figure 6 . -Pressure distributions on a parabolic-arc 
body of revolution of fineness ratio G at subcritical Mach 
numbers as indicated by present theory. 


theory are shown in figure 7. These results are 
repeated in figure 8 together with the correspond- 
ing theoretical results indicated by the second- 
order approximation given by equation (32), and 
the experimental results of Drougge from refer- 
ence 14. The latter results were obtained from 
tests with a model that consists of the forward 5/6 
of a parabolic-arc body of revolution of fineness 
ratio 6. It was supported in the wind tunnel by 
a circular cylinder of the same diameter as the 
base of the body. Although no corrections have 
been applied for the effects of the abbreviated 
length of the body, it is clear that the values for 
Cp on a complete body should be somewhat more 
negative than those on the partial body tested by 
Drougge. Theoretical considerations indicate that 
the magnitude of this effect is about 0.007 at the 



Figure 7.- -Pressure distributions on a parabolic-arc 
body of revolution of fineness ratio 6 at subcritical Mach 
numbers as indicated b} T linear theory. 
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middle of the body, and about 0.001 at the nose 
of the body. A more detailed investigation of this 
matter can be found in reference 10 by Van Dyke. 

Two conclusions emerge from a comparison 
of the theoretical and experimental results illus- 
trated in figure S. They are that slender-body 
theory is capable of yielding reliable results for 
subsonic flow past a parabolic-arc body of revo- 
lution of fineness ratio 6, and that the simple 
result provided by linearized compressible flow 
theory remains a useful approximation for all 
Mach numbers up to the immediate vicinity 
of the critical 

Theoretical pressure distributions for a para- 
bolic-arc body of revolution of fineness ratio 6 
could also have been calculated, and with less 
effort, by application of the similarity rule given 
by equations (9) and (12) to the results already 
calculated for the parabolic-arc body of fineness 
ratio 10. As noted in the course of the preceding 
discussion of the similarity rule, the results 
obtained by application of equations (9) and (12) 
do not agree perfectly with those obtained by 
direct calculations. It is the purpose of the 
present discussion to illustrate the magnitude 
of these differences by an examination of the 
specific case mentioned above. Figure 9 shows 



Futtrk 9.- — Variation of pressure coefficient with Much 
cumber at various stations on a parabolic-arc body of 
fineness ratio G as obtained by direct calculations using 
the present theory, by use of the similarity rule, and by 
experiment. 


the variation of C p with frec-stream Mach number 
at several stations along the length of the body 
as indicated by the direct calculations for the 
body of fineness ratio 6, by application of the 
similarity rules together with the results calcu- 
lated for the body of fineness ratio 10, and by the 
experimental results of Drouggc (ref. 14) for a 
body of fineness ratio 6. The stations for which 
results are included on this graph are selected 
so as to coincide with the location of the orifices 
in the model tested by Drougge. It can be seen 
that the results obtained by application of the 
similarity rule are very close to those obtained 
by direct calculation over substantial portions of 
the body, but that important differences occur in 
the vicinity of the middle of the body as the 
frec-stream Mach number approaches the critical 
and the local velocities approach the speed of 
sound. Perhaps the most striking difference 
concerns the frec-stream Mach numbers at which 
the two sets of theoretical results terminate. 
In the direct calculation, the results terminate 
with the attainment of the critical Mach number, 
whereas the results obtained by application of 
the similarity rule terminate at Mach numbers 
well below the critical. If it is presumed that 
results appropriate for supercritical Mach num- 
bers would contain shock waves, and that the 
bodies would experience a pressure drag, it is 
clear that the results obtained by application of 
the similarity rules would be in serious error 
since, in the present example, such effects would 
begin to appear at a Mach number of about 
0.845 rather than about 0.88. It is evident, 
furthermore, that the magnitude of these effects 
would be less if the fineness ratios of the two 
bodies were more nearly the same, but greater 
if the fineness ratios were more different, and of 
opposite sign if the similarity rule were applied 
to calculate the pressure distribution on a thin 
body using known results for a thicker body. 

SUPERSONIC FLOWS 

APPROXIMATE SOLUTION OF EQUATIONS 

The procedure described in the preceding sec- 
tion will now be applied to the analysis of super- 
sonic flows. Thus, introduce the symbol \ H as 
an abbreviation for the coefficient of p xx 

X*=A/.*- 1 +MJ 1 +hu > 0 

* ' CO 

( 34 ) 
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and rewrite equation (4) as follows: 

^H<Pxx — ~ <Pr — <Prr— 0 (35) 


If it is again assumed that is neither zero nor 
infinite and that it varies sufficiently slowly that 
its derivatives can be disregarded, the problem 
is equivalent to that encountered in linearized 
supersonic slcnder-body theory. The solution 
for <p for points on the surface of a smooth slender 
pointed body of revolution is given by (see, 
e.g., ref. 8, p. 190) 


— u 


-VV S'(£)d£ 


VO ** 1 


-£) 2 -x„r* 


(36) 


where the subscript II denotes that the solution 
refers to equation (35) of hyperbolic type with 
X H constant. The corresponding expression for 
u H can be approximated for points on the surface 
of a smooth slender body by 


U m 


S"(x) , V? 
4tt 4ttx 2 


2irJo X—i 


dt 


S”(x) 

4 x 


In Xjy+/# (a?) 


(37) 


where J H (x) can be interpreted as the expression 
for u/Un indicated by linear theory for a free- 
stream Mach number of \2. Differentiation 
yields 


diUiifUJ S'" (x),^ , df H (x) 
dx ~ 4tt ln X//+ dx 


(38) 


constant for flows that are supersonic everywhere 
is to use the value for u at that is provided 
by the solution of the equations of transonic flow 
theory for flow around a cone tangent to the nose 
of the body. Practical difficulties arise in the 
application of this procedure, however, because 
first, the solution for the cone cannot be written 
in a simple form, and second, the numerical tables 
and charts in which the results of numerous specific 
cases are summarized are generally inadequate for 
the present purposes because of the lack of suffi- 
cient calculated results for slender bodies for low 
supersonic Mach numbers. These comments 
apply not only to the fairly recent results of 
Oswatitscli and Sjodin (ref. 15) and Shen (ref. 16) 
based on the simplified equations of transonic flow 
theory, but also to the older and more extensive 
results based on the complete equations of com- 
pressible flow theory and presented in references 
17 and 18. An alternative method that avoids 
this difficulty and parallels the procedures em- 
ployed for the approximate solutions of other 
problems in this paper and in reference 4 is to 
evaluate the constant by assuming that u equals 
u H at the point on the forebody where S n {x) 
vanishes, since u H is independent of \ H and hence 
of Mach number at this point; thus 

where S n (x)=0 (40) 

C CO 

The calculation of u on the surface of a body of 
revolution of specified shape in supersonic flow 
is thus reduced to the solution of an ordinary 
differential equation of the form 


If, in the same manner as described for the sub- 
sonic case, M w 2 - 1 -\-Jcu is restored in place of 
X/, so that, in effect, the local value for X// is 
used at each point and the subscript II on u is 
dropped, equation (38) becomes 


dx 


S"’(x ) 
4t 


ln(MJ-l+Jcu) 


(Ifni* ) 

dx ' 


(39) 


As in the previous discussion of subsonic flows, 
the resulting relation for u on the body surface is 
a nonlinear ordinary differential equation of the 
first order, and it is necessary to supply a value 
for u at some point along the body in order to 
evaluate the constant of integration. Perhaps 
the most logical method for the evaluation of this 


fx =F ^ ^ 

subject to the additional condition given by equa- 
tion (40). This problem is identical in form to 
that encountered at the corresponding point in 
the discussion of subsonic flows, and can be solved 
numerically in a similar manner. Once the values 
for u on the surface of the body are determined, 
the corresponding values of C p follow directly 
upon insertion of the results into equation (3). 

APPLICATION TO PARABOLIC-ARC BODY 

As an application of the foregoing theory, con- 
sider supersonic flow past a parabolic-arc body of 
revolution of length l and maximum diameter d. 
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The expressions for the ordinates 7? , the cross- 
section area S , and the location of the points at 
which S"(x) vanishes are given by equations (25), 
(26), and (27). Substitution of these relations 
into equations (37) and (39) leads directly to the 
following differential equation for u on the surface 
of a parabolic-arc body of revolution in super- 
sonic flow: 


latter results are calculated by use of the following 
expression derived from equations (3), (26), and 
( 37 ) with X H replaced by MJ— 1 . 

+24t 2 ^(2-3 y)-/?' ! (44) 


\n(MJ-X+ku) 

+s{If 1 "[' O-O]- 12 '’ K 2-3 0} (42i 

The constant of integration is evaluated by im- 
posing the condition derived from equations (26), 
(27), and (40) that 

r f=-2 v '3r 2 at *=4(l~y) (43) 

The problem described by equations (26), (42), 
and ( 43 ) has been solved numerically, by applica- 
tion of Milne’s method described in the preceding 
section on subsonic flows, for several supersonic 
Mach numbers for a body of fineness ratio l/d of 
10. The resulting values for u, A/29, and C P are 
tabulated in table III. The results for C p are 
also shown graphically in figure 10. The corre- 
sponding results indicated by linearized compress- 
siblc flow theory are shown in figure 11 * The 


The results presented in table III and figures 
10 and 11 are repeated in figure 12 together with 
the corresponding results indicated by the follow- 
ing expression given by Van Dyke in reference 10 
for the second-order approximation to the solu- 
tion of the equations of transonic flow theory for 
slender parabolic-arc bodies of revolution: 


c -c r ,r 3 nfg'Y.g/gyi 

+ MJ - 1 l 6 2 \ T / ' 4 \ T / _ 


As noted previously in the discussion of subsonic 
flows, the quantity 7+1 that appears in the result 
given by Van Dyke has been replaced by 
M w 2 ( 7 + 1 ) to conform with the present formula- 
tion of the equations of transonic flow theory. 

It can be seen upon examination of figures 10 , 
11 , and 12 that the results obtained by applica- 
tion of the present theory arc again qualitatively 
similar to those indicated by first- and second- 
order theories and that the main differences occur 



Figure 10. — Pressure distributions on a parabolic-arc Figure 11— Pressure distributions on a parabolic-arc 

body of revolution of fineness ratio 10 at Mach numbers body of revolution of fineness ratio 10 at Mach numbers 

above the upper critical as indicated by present theory. above the upper critical as indicated b^ linear theory. 
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(a) Af eo = I.2I (b) M m = 1.22 (c) M ro = 1.23 

(d) = 1.25 (e) Mo, = 1.30 W M ro = 1.40 

Figure 12. — Pressure distributions on a parabolic-arc body of revolution of fineness ratio 10 at Mach numbers above the 
upper critical as indicated by present theory and by other theoretical results. 

501377 — 59—3 
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in the regions where the local Mach number ap- 
proaches unity. The most prominent difference 
occurs near the tail of the bodies where the first- 
and second-order theories indicate that the flow is 
subsonic. 

This result is contradictory to the basic assump- 
tion that the flow is supersonic, but it is never- 
theless qualitatively correct, since more exact 
investigations also indicate that the flow is sub- 
sonic at the tail of a body of revolution that termi- 
nates like a cone. There is no assurance that the 
smooth deceleration through sonic velocity indi- 
cated by the first- and second-order theories is 
correct, however, and further investigation appears 
necessary for the clarification of this detail. On 
the other hand, the results indicated by the present 
theory terminate when u reaches the value 
associated with the occurrence of sonic 
velocity. When this situation occurs, the inequal- 
ity specified by equation (34) is no longer satisfied, 
the argument of the logarithmic term of equation 
(39) vanishes, and dujdx increases logarithmically 
to infinity. If the solution for the remainder of 
the body is required, it is necessary to consider the 
nature of the shock wave that must be present in 
the vicinity of the rear of the body and to use 
equations that are not restricted in application to 
supersonic flows. As a result, further discussion of 
this matter is deferred to the subsequent section 
on flows with free-stream Mach number near 1. 
In any case, the entire question is of rather an 
academic nature for parabolic-arc bodies in the 
Mach number range considered in the preceding 
discussion because (a) the region involved is only 
a small fraction of the body length, (b) significant 
viscous effects are generally present in this region, 
and (c) the information is frequently not needed 
because the bodies used in practice or investigated 
in the wind tunnel are cut off forward of this 
region. 

Xo experimental pressure distributions for a 
parabolic-arc body of revolution of fineness ratio 
10 are available for comparison with the theoreti- 
cal results shown in figure 12. The maximum 
test Mach number of the investigation reported 
in reference 11 is 1.20, whereas the lowest Mach 
number for which the present theoretical results 
indicate supersonic flow at the nose is 1.21, The 
comparison with first- and second-order theory 
should suffice, however, for an estimate of the ac- 
curacy of the present results for supersonic flows. 


FLOWS WITH FREE-STREAM MACH NUMBER 
NEAR 1 

APPROXIMATE SOLUTION OF EQUATIONS 

The analyses of subsonic and supersonic flows 
given in the preceding sections have started by 
introduction of a symbol X for the coefficient of 
(p xx and the assumptions that X is neither zero nor 
infinite and that it varies sufficiently slowly that 
it can be regarded as a constant in the early stages 
of the analysis. Since the results so obtained 
terminate if X=0, or physically if a point is ap- 
proached at which the velocity is sonic, it is evi- 
dent that additional considerations are necessary 
to permit the study of flows with free-stream Mach 
number near 1 in which the acceleration from sub- 
sonic to supersonic velocities is an essential feature. 

The technique adopted is to introduce the sym- 
bol X/> as an abbreviation for the coefficient of 
<p z rather than <p XX) thus 

A,. MJ g>0 (46) 

whence equation (4) may be written as follows: 

~ tpr-btPrr Xp<£ r = (il/^ 2 1 )<Prz = fp (d7) 

It is assumed, once again, that \ P is nonsingular 
and that it varies sufficiently slowly that it can 
be considered as a constant in the initial stages of 
the analysis. The resulting differential equation 
is linear and is of elliptic or hyperbolic type de- 
pending on whether the free-stream Mach num- 
ber, rather than the local Mach number, is less 
than or greater than unity. At free-stream Mach 
number 1, the term f P vanishes and equation (47) 
reduces to a partial differential equation of para- 
bolic type that is familiar from the study of heat 
conduction. This appears to be reasonable 
inasmuch as properties of equations of parabolic 
type are intermediate, in a certain sense, between 
those of equations of elliptic and hyperbolic type. 
At free-stream Mach numbers slightly different 
from 1, the term f P remains, but is small. The 
present analysis of such flows is therefore also 
based on considerations that are normally applied 
to equations of parabolic type. 

The starting point is the following relation for 
<p P derived from equation (47) with X P positive 
and constant, the boundary conditions stated in 
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equation (1), and the form of Green’s theorem 
associated with the left side of equation (47): 


4 TTjo X i 


1 f*2ir ('<» ri 

~ de pdp <r P f r tla (48) 

h P J o Jo J- co 


where 


(T P = 


x P 

"4t r(x—& 


— 2r/j cos (0— 0)1 


The corresponding expression for u P follows by 
differentiation. It can be approximated for points 
on the surface of a smooth slender body by 

Ur S"{x) 1 1 r S"(£)-S"{J) 

U m 47 r 47r./“ 47 tJ 0 r — £ s 

«•> 

where 


C r =Euler’s constant — 0.577215GG5 . . . 


If the free-stream Mach number is unity, f P 
vanishes and u p can be calculated directly. The 
resulting expression is precisely that proposed by 
Oswatitsch and Keune in references 19 and 20 for 
the calculation of u on the surface of the forepart 
of an arbitrary slender, pointed body of revolution 
at free-stream Mach number 1. They have shown 
that the results arc in remarkably good agreement 
with those measured on the front half of a para- 
bolic-arc body having a thickness ratio t of %. 
In this comparison, the parameter X P was arbitrar- 
ily equated to 


Xp a 2 (7+1) lj V a 2 2 

in the notation of the present paper. In the in- 
vestigations mentioned above as well as further 
studies reported in references 21 through 26, the 
parameter \ P is regarded throughout as a con- 
stant, and various means are proposed for the 
selection of an appropriate value. It appears, in 
general, that the results obtained in this way are 
remarkably accurate if the resulting values for 
bu/bx are indeed reasonably constant. If, on the 
other hand, bu/bx varies sufficiently along the 


body, no choice of a single value for \ P will lead 
to a useful result. Tins point, which has already 
been discussed in some detail for two-dimensional 
flows in reference 5, is developed further in the 
course of the following discussions. Some criti- 
cisms of the above procedure, although principally 
from a different point of view, have appeared in a 
note by Miles (ref. 27). 

If the free-stream Mach number is not unity, 
equation (49) is an integral equation, and it might 
appear that little progress toward a solution has 
been made. If attention is confined to Mach 
numbers near unity, however, it is only necessary 
to approximate well locally and it is sufficient 
to substitute X P /k for <px or bu/bi in the triple 
integral. The integral can then be evaluated and 
the following relationship results: 


u P \-\fJ 

U m MJ{ 7+1) 


S''(x) \ P St r 
4tt 11 4ttj 


MttJo i-f 




(51) 


If, once again, k (bu/bx) is restored in place of 
A P , the subscript P is dropped, and use is made of 
the following relation between bu/bx and du/dx 
along the surface of the body 


d / u 
dx \U c 


gte > 


dJR 
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srsr 

4irS 


(52) 


a nonlinear ordinary differential equation is 
obtained for u oil the body surface. It is 
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This equation is of the first order and the same 
general form 

( 54 ) 


as encountered in the analysis of subsonic and 
supersonic flows in the preceding sections, and it 
is again necessary to specify some condition for 
the evaluation of the constant of integration. 
Two classes of cases are to be distinguished. One 
class includes the foreparts of bodies with a convex 
corner, such as a cone-cylinder, for which the 
necessary condition is supplied by the fact that 
(he velocity must be sonic, that is, 


u 1-3/, 2 

U m MJ(y+ 1) 


(55) 


follows that all of the coefficients of a Taylor 
series expansion for u at the point where S"(x) is 
zero can be determined, the only requirement 
being that, at each step, the derivative of next 
higher order is finite. The Taylor scries suffices 
for the determination of the solution of equation 
(53) in the neighborhood of the point where S"(x) 
vanishes. The remainder of the solution can be 
determined by application of standard methods 
such as that of Milne described in the preceding 
section on subsonic flows. 

APPLICATION TO CONE-CYLINDER 

A simple application of the foregoing theory is 
provided by consideration of flow with frce-slroam 
Mach number near 1 past a slender circular cone- 
cylinder of maximum diameter d and cone length 


at the shoulder. The other includes smooth 
bodies along which the velocity accelerates 
smoothly through sonic velocity and for which 
there is no point for which a value for u can be 
specified by a priori considerations. The pro- 
cedures followed in the preceding sections and in 
reference 5 suggest that u be equated to u P at the 
point on the forepart of the body where S rr { x) is 
zero, since, for Mach numbers near unity, u P is 
independent of at this point. This point lias 
the additional property, evident from equations 
(53) and (55), of corresponding to the station 
toward which the sonic point approaches as the 
thickness ratio of the body becomes vanishingly 
small. 

This procedure does not suffice to determine a 
solution in the present case, however, because 
equation (53) is singular at this point and an 
infinity of integral curves pass through the value 
v = u F at the point where S /; (x) vanishes. Of all 
these curves, however, only one is analytic (all 
derivatives finite) at this point, and selection of it 
suffices to determine a unique solution. That 
this is so can be seen by consideration of the series 
of expressions obtained by successive differentia- 
tion of equation (53a) with respect to x. After 
an arbitrary number n of differentiations, an 
expression is obtained for d n ufdx n that involves 
derivatives of u with respect to x up to the order 
n-f 1. The tt + 1 derivative is, of course, unknown, 
hut it is multiplied by S”(x) and the product 
vanishes at the point where S"(x) is zero, pro- 
vided always that d n+l ujdx n+l is finite. It thus 



1 1 2 as illustrated in figure 13. The ordinates R 
of this body are given by 


R~d j~ Tjr for 0<+<7/2 
R=^=t ^ for xd>l / 2 


(56) 


The variation of cross-section area S with x along 
the conical part of the body is thus 


S= J !rR 2 — 7T T 2 X 2 


(57) 


Details of the solution. — Substitution of the 
expression for S(x) given in equation (57) into 
equation (53b) leads to the following differential 
equation for u along the surface of the cone: 


d_V u V+-1 
ds\jU + MJ( 7+1) 
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Introduction of the new variable G(x) defined b_y 
MJ-l 
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(59) 
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leads to an equivalent differential equation 
dO T 4r 1 2 G 

S = Lv.«<* + l)rV’J “ P Y <60) 

which can be solved by separation of variables. 
The result is 

G= ~2 ln \_ Cl ~MJ{y+ 1) r i e c '] (G1) 

The corresponding expression for u along the 
surface of the cone follows from equation (59) 
and is 

- r H c '-M7£ W] ™ 

The constant of integration C 2 that appears in 
equations (61) and (62) is evaluated by use of the 
condition mentioned in the preceding section that 
the local velocity must be sonic at the shoulder, 
that is, the relation given by equation (55) holds 
at x~lj 2. This procedure leads to the following 
result: 

<*-G) } + C 2 (7+l)d< c ] (63) 


Substitution of this result for C r in equation (62) 
and a slight rearrangement of terms yields the 
following expression: 



The corresponding expression for the pressure 
coefficient C p on the surface of the cone follows 
from equation (3) and is 


P ~Mj(y+ 1 ) 


— 2r 2 In 


TX 

m 



A plot of the values for C v that result upon appli- 
cation of equation (65) to the special case of a 
cone-cylinder of semiapex angle t of 1/10 in a flow 
with frec-stream Mach number 1 is shown in 
figure 14. 
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Figure 14.— Pressure distribution on a cone-cylindcr at 
Mach number 1, as indicated by equation (G5). 

Examination of the result given b}- equation (65) 
reveals that the similarity rule for the surface 
pressures given by equations (9) and (12) is not 
satisfied exactly. It is evident, however, that 
perfect agreement would be attained if the follow- 
ing approximation were introduced: 


4-G^yi 

d f w 2 (7 + 1 ) T 2 e c M m 2 (7 + 1 ) r 2 c c 


(06) 


It follows, therefore, that the similarity rule is 
very nearly satisfied over most of the length of 
the cone, and that the anticipated failure of the 
similarity rule in the vicinity of the shoulder is 
effectively confined to a very small region. 

Further examination of the result given by 
equation (64) reveals that the condition that 
dufbx is greater than zero specified by equation 
(46) at the outset of the analysis is not satisfied at 
the apex. This is evident from the following ex- 
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prcssion derived from equations (52), (57), and 
(64) 

r & /v_\ i s>s " 

Ld-r V '»/ Jk dx \U a ) R 4 tS 


4t ’G k ) 

MJ(y+ l)rV+4[.-(i)’]} 


from which it can be seen that c Xw/Su? vanishes at 
This result suggests that the range of use- 
fulness of equation (65) has probably been ex- 
ceeded in the vicinity of the apex, and recalls an 
example in two-dimensional flow presented in 
reference 4 in which bujdx vanishes and then 
becomes negative along the rear part of a par- 
ticular airfoil. 

Determination of correction for the vicinity of 
the nose. — The calculation of the pressure distri- 
bution at a Mach number near 1 on a body of 
such a shape that du/da? is positive along part of 
body and negative along the remainder cannot be 
accomplished by direct application of any one of 
the relations developed in the preceding sections. 
On the one hand, the considerations based on a 
differential equation of parabolic type that lead 
to equation (53) permit the analysis of flows that 
pass through sonic velocity, but fail when d-u/dx 
is zero. On the other hand, the considerations 
based on a differential equation of elliptic or hy- 
perbolic type that lead to equations (20) and (39) 
permit the analysis of flows in which c)u/ c);r passes 
through zero even though the free-stream Mach 
number may be unity, but fail if the local velocity 
is sonic. The breakdown is associated in each 
case with the fact that the basic partial differential 
equation, that is equation (16), (35), or (47), as- 
sumes a degenerate form when X is zero. The 
pressure distribution along a body having the 
properties described above can, nevertheless, be 
calculated by considering the solution in sections 
and joining together the various results in such a 
way that the failings associated with vanishing X 
are avoided. This procedure, which is much 
simpler than a complete re-analysis of the problem 
from a sufficiently general point of view to en- 
compass the entire problem in a single sweep, is 
exactly the same as that already employed in 
reference 4 in the analysis of the related problem 
in two-dimensional flow mentioned above. 


The flow in the vicinity of the apex of a slender 
cone-cylinder is subsonic in nature, hence the 
further investigation of the pressure distribution 
on this part of the body is based essentially on 
equation (20). It is necessary, however, to add 
some qualifying remarks about the expression to 
be used for u u the solution for zero Mach number, 
since the discontinuity in the slope of the surface 
at the shoulder violates the conditions imposed in 
the derivation that the body be smooth. The 
particular case of a slender cone-cylinder in in- 
comprehensible flow, or in linearized subsonic flow, 
lias been investigated by several authors and ap- 
proximate solutions involving various degrees of 
refinement and complexity are presently available. 
It is not necessary to employ the more refined 
approximations in the present application, how- 
ever, since only the part of the result that pertains 
to the forward part of the cone is to be used. A 
simple expression that is adequate for the present 
purposes is one that follows from equation (17) 
upon carrying out the indicated operations, ex- 
panding the terms containing t in a power series, 
and disregarding terms proportional to the third 
or higher power of r. The resulting expression 
has been given by Lai tone in reference 28, and 
can also be obtained directly from somewhat 
more general expressions given in references 29 
and 30, or on page 311 of reference 6. It is 



The corresponding expression for the pressure 
coefficient C P can he determined by substitution 
of this relation into equation (3). A measure of 
the accuracy of this approximation can he had by 
comparison with the results of a much more 
refined analysis of incompressible flow past a 
slender cone-cylinder given by Fraenkcl in refer- 
ence 31. The results of such a comparison are 
illustrated in figure 15 for the special case of a 
cone-cylinder having a semiapex angle r of 1/10. 
It can be seen that, the two theories are in very 
good agreement over the forward portion of the 
cone. 

Differentiation of equation (68) in accordance 
with the procedure associated with equations (18) 
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Figure J5. — Pressure distribution on a cone-cylinder in 
incompressible How. 


and (19) results in the complete loss of X* in the 
analysis. The subsequent integration is immedi- 
ately possible, therefore, and the resulting ex- 
pression differs from equation (68) only in that 
the term involving is replaced by an arbitrary 
constant of integration, thus 





0 ~m) 


+C t (69) 


The value for the constant of integration C t is to 
be selected so that the value for u given by equa- 
tion (69) is equal to that given by equation (64) 
at the point where the two results are joined. 
The latter point is determined by selecting the 
point at which the two values for du/dx obtained 
by differentiating equations (64) and (69) are 
equal. It is found, upon following these pro- 
cedures and using the approximate relation stated 
in equation (66), that the point at which the two 
results are to be joined is situated at 


jr 1 

W~ 3 


(70) 


and that the expression for the constant of inte- 
gration is given by 

r 3 , i 2 , r 4 i 

MJ (7+1) 4 7 2 1,1 Uc*'' '(7+1 )" tVU 

(71) 


These results lead to the following expression for 
C p upon assembly and insertion into equation (3) : 


+ r2 '" {[.T/„ 2 (7+1) rv] [(772) (' “7/2)]} 
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IP J 


-r 2 (72) 


It follows furthermore that equation (72) is to be 
used to calculate the pressure distribution for the 
forward third of the cone, and equation (65) is to 
be used for the remainder of the length of the cone. 

The magnitude of the difference between the 
two results and the manner in which they join 
together are illustrated in figure 1C for the special 
case of a cone-cylinder of semiapex angle r of 1/10 
in a flow with free-stream Mach number 1. This 
figure shows the values for C v indicated by equa- 
tion (72) for the front portion of the cone together 
with those indicated by equation f65) for the en- 
tire length of the cone. It is evident from this 
comparison that the two results can be joined 
together in a manner that is acceptable for most 
practical purposes. Although the difference be- 
tween the two results is not large over most of 
the forward third of the cone, the result indicated 
b} r equation (72) confirms the general motivation 
for the determination of a correction for this 
region since the values for bu/bx are negative in 
the vicinity of the apex. It may be noted, fur- 
thermore, that the results indicated by equation 
(72) are in perfect agreement with the similarity 
rule indicated by equations (9) and (12). This is 
a consequence of the systematic use of the approx- 
imate relation given in equation (66). 

Comparison with experimental and other theo- 
retical results. “Experimental results for C p on 
the surface of a cone-cylinder having a semiapox 
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Figure 16. — Illustration of nose correction to pressure 
distribution on a cone-cylinder at Mach number 1. 


angle of approximately 7° (r — 0.1225) have been 
given for a range of Mach numbers around unity 
by Page in reference 32. Since it was considered 
particularly important in Page's investigation to 
both minimize and ascertain the degree to which 
the results arc influenced by wind-tunnel wall in- 
terference, the tests were conducted in two wind 
tunnels of different size, one having a test section 
14 feet square, the other 2 feet square, using a 
small model, the conical part of which was 5.50 
inches in length. The results for freo-stream Mach 
number 1 are presented in figure 17. The differ- 
ence between the two sets of experimental results 
is attributed by Page to wind-tunnel wall inter- 
ference. Page goes on to show that substantial 
interference effects remain even in the data ob- 
tained in the larger wind tunnel, and finally de- 
termines, by both theoretical and experimental 
considerations, a curve for the pressure distribu- 
tion corrected for interference. The latter results 
are indicated in figure 17 by the flagged symbols. 
The corresponding theoretical pressure distribu- 
tion calculated by use of equations (65) and (72) 
is indicated by a solid line. Also included in this 
figure is a plot of the pressure distribution indi- 
cated by the numerical approximate solution of 



Figure 17. — Pressure distribution on a cone-cylinder at 
Mach number 1, as indicated by present theory, by 
other theoretical results, and hy experiment. 


Yoshihara (ref. 33) 1 for a cone-cylinder having a 
semiapex angle of 1/10 radian transformed so as 
to bo appropriate for a 7° cone by application of 
the similarity rule given by equation (12). It can 
be seen that the pressure distribution given by 
the present theory is distinctly different from that 
given by Yoshihara, but is similar in form to the 
experimental pressure distributions, although dis- 
placed somewhat. The displacement is consider- 
ably less for the tests conducted in the large wind 
tunnel than in the small wind tunnel, and nearly 
disappears when the results of the tests in the 
large wind tunnel are corrected for interference. 

The foregoing results may be compared with 
those obtained by application of the simple line- 
arized theory of Oswatitsch and Keune (refs. 19 
and 20) for sonic flow past a slender body of revo- 
lution. In this theory, X P is replaced by a con- 
stant, such as that indicated by equation (50), 
and the corresponding expression for the pressure 

i A correction has been applied to Yoshihara’s results to allow for a sign 
error in the quadratic term of the expression for C v . 
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Figure 18. Comparison of results indicated by equation 
(73) with those indicated by the present theory for tin- 
pressure distribution on a cone-cylinder at Mach 
number 1 . 


distribution on the conical part of a slender cone- 
cylinder follows directly from equations (51) and 
(3). It is 


r,— 


(73) 


if equation (50) is used to evaluate \ P . Figure 18 
shows a plot of the numerical results for a cono- 
cylinder having a semiapex angle of 7° together 
with the present theoretical results and the experi- 
mental results of Page, as corrected for wall inter- 
ference. It can he seen that the results indicated 
by equation (73) arc of more or less the correct 
order of magnitude, hut fail to predict, even quali- 
tatively, the proper behavior in the vicinity of the 
shoulder. It should he observed, furthermore, 
that no essential improvement would result if 
were replaced by some other constant that than 
indicated by equation (50) since such a change 
leads only to a uniform change in C p along the 
entire length of t lie cone. This difficulty of the 
linearized theory for sonic (low is to he expected 


because the region of dependence for the solution 
at an arbitrary point is confined to the part of 
space situated upstream of this point. As a re- 
sult, the influence of the part of the cone-cylinder 
situated between the arbitrary point and the 
shoulder is not included in the solution, although 
it should he since the flow in this region is sub- 
sonic. The present theory, on the other hand, 
includes the influence of this region in the final 
result since the constant of integration associated 
with the differential equation that results upon 
replacement of X/. by k(du[df) is specified at the 
shoulder and the integration proceeds in the up- 
stream direction. Tn this way, the approximate' 
expression for C p indicated by equation (65), and 
therefore also that indicated by equation (72), de- 
pends on conditions along the entire conical part 
of the cone-cylinder, and the influence of regions 
situated both upstream and downstream from an 
arbitrary point are included in the result. 

Theoretical results for the pressure distribution 
on a slender cone-cylinder in flows with high sub- 
sonic or sonic free-stream velocity have also been 
given recently by Tvusukawa in reference 34. 
These results were obtained by application of an 
approximate relation between axisymmetric and 
two-dimensional flows given previously by the. 
same author in reference 35. In this way, the. 
solution for a cone-cylinder is related lo that for 
an airfoil having ordinates proportional to r 2 fol- 
lowed by a straight section extending downstream 
to infinity. Tvusukawa obtains an approximate 
solution for the latter problem by considering a 
flow that separates from the shoulder of the airfoil 
and calculates the pressure distribution on the air- 
foil surface with the aid of a pressure correction 
formula. This formula, which was given in a 
third paper by the same author (ref. 36), relates 
the pressure at a given point in a compressible 
flow to that at the same point in an incompressible 
flow past the same airfoil. This procedure has 
been applied to calculate the pressure distribution 
on a cone-cylinder having a semiapex angle of 7° 
and the results are presented in figure 19 together 
with the present theoretical results and the experi- 
mental results of Page. The latter results are 
again corrected for wall interference, as noted in 
the discussion associated with figure 17. It 
appears from this comparison that the values for 
C p calculated by (lie method of Tvusukawa are too 
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Fkjcke U). -Comparison of results indicated by Kusu- 
kawa (ref. 34) with those indicated by the present theory 
for the pressure disl ribution on a cone-cylinder at Mach 
number 1 . 


large at all points along the length of (he cone. 
Of the three principal steps in the derivation of 
Kusukawa’s results, the only one that has not 
been checked by a number of comparisons with 
experimental or other theoretical results, and is 
thus the most likely source of the discrepancy, is 
the approximate relation between axisymmetric 
and two-dimensional flows. It is not likelv that a 
significant part of the discrepancy is a conse- 
quence of the use of the assumption that the flow 
separates from the shoulder instead of performing 
a Prandtl-Meyer expansion around the corner, 
since it has been shown in references 3G, 37, and 5 
for the case of a wedge profile that the nature of 
the flow behind the shoulder has very little effect 
on the pressure distribution on the part of the air- 
foil ahead of the shoulder. Neither is it likely 
that a significant part of the discrepancy can be 
traced to the use of the pressure-correction 
formula. This formula, which was derived by 
Kusukawa by application of the WKB method of 
approximation to the equations of transonic flow 


theory expressed in terms of hodograph variables, 
lias also been found independently (except for 
replacement of y+1 by A/ co 2 ( 7 +l) as a result of a 
corresponding change in the fundamental equa- 
tions of transonic flow theory) by application of 
much different considerations and is presented as 
equation (28) in reference 4. This result is 
evaluated in references 4, f>, and 36 by numerous 
comparisons with other theoretical and experi- 
mental results. Of the various examples, the one 
that is most closely associated with the present 
discussion is a comparison given in reference 5 
with the results indicated by the exact solution of 
the equations of transonic flow theory for sepa- 
rated flow past a wedge profile published recently 
by ITelliwell and Mackie in reference 37. This 
comparison shows that the numerical values for 
the pressure distribution on the surface of the 
wedge are hardly distinguishable, whether calcu- 
lated by use of the exact or approximate formulas. 

Some insight into the range of Mach numbers 
surrounding unity for which the present theory 
might he expected to give useful results can be 
gained by noting that the indicated variation of 
C p with M w is essentially proportional to M n - 1 
and that the constant of proportionality is very 
nearly equal to that- associated with the known 
invariance of local Mach number A/ with changes 
in free-stream Mach number at A/« = l, that 
is, with the relation 



The exact relation for dCJdM « for A/„ — 1 for 
isen tropic flow that corresponds to equation (74) is 



2 

7+1 




(75) 


The approximate relation indicated by the present 
theory is 



(76) 


except at the shoulder where it is 

(dC p \ = _4 
7+1 


(77) 


All three relations are nearly equivalent, since 
2Cy( T +l) and 2r 2 are very small in comparison 
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with 4 /( 7 + 1) for slender bodies. It follows, 
therefore, that the variation of C p witli indi- 
cated by the present theory is sufficiently accurate 
for all engineering purposes within at least a small 
range of Mach numbers surrounding A/ w = 1. It 
should be noted that the range of Mach numbers 
for which the above results may be expected to 
apply is much smaller than in two-dimensional 
flows, since Guderley has indicated (ref. .38, p. 29G) 
that the leading term in the expansion of the 
deviation of M from its value for A/«, = l is pro- 
portional to (A/„ — l) 3 for planar flows and 
(A/* — 1) 5/3 for a Asymmetric flows. 

Pressure drag. -The pressure drag of a slender 
cone-cylinder at free-stream Mach number 1 and 
vicinity can be found by integration of the general 
relation given by equation (13) with the expres- 
sions given by equations (57), (05), and (72) 
substituted for S and C p . This procedure leads to 
the following expression for the drag coefficient 
referred to maximum cross-section area if equa- 
tion (G5) is used for C p along the entire length of 
the cone: 

C — - - ^ —2 -y. a ~* 

(yfV) kr 2 (//2) 2 l M “ 2(7+,) 

+ r 2 ^-1-f ln Vf j {y+ i) T vcJ (' 8 ) 


If equation (72) is used instead of equation (G5) 
to calculate C p along the forward third of the cone, 
the value for C D is smaller by an amount given by 




(79) 


so that 


M 2 — 1 


7 1 1 ) 


,r_i , , 04 

L 2 ^ 27MJ(y\\)rY_ 


(80) 


Comparison with equations (9) and (14) reveals 
that both of these results for C D satisfy the sim- 
ilarity rule exactly. 

The variation of C D with 7 for- free-stream Mach 
number 1 indicated by equation (80), which 
reduces to the following simple form for ^f ao ~ 1 
and y — 1 .4, 

(81) 



Fiourk 20. -Viiriation of pressure-drag coefficient witli 
thickness ratio for a cone-cylinder at Mach number I, 
as indicated by present theory, by oilier theoretical re- 
sults, and by experiment. 

is shown graphically in figure 20. Also included 
in this figure are the corresponding theoretical 
results indicated by the following equation deter- 
mined by combination of the similarity rule and 
Yoshihara’s result for a cone-cylinder having a 
semiapex angle of 1/10 radian: 

^,= -7*0.55-1 4 1 nr) (82) 

and by the following equation given by Miles in 
reference 39. 

fh=-r 2 (0.091 -f 4 In 7 ) (S3) 

The latter formula is derived by an approximate 
procedure based on matching a solution of La- 
place’s equation in t wo dimensions for points near 
the axis with the asymptotic solution of Guderley 
and Yoshihara (refs. 40 and 41) for points at great 
distances from the axis. The solutions are joined 
at the sonic line, and the drag is determined by 
application of momentum methods without an 
intermediate calculation of the pressure distrihu- 
( ion . Three data points for a cone-cylinder having 
a semiapex angle of 7° are also included in figure 
20. They represent the values obtained from the 
experimental investigation of Page reported in 
reference 32. In accordance with the notation 
already employed in the display of Pago’s results 
in figure 17, the square symbol represents the 
result obtained from tests in a 2-foot wind tunnel, 
the plain circle represents the result obtained in 
a 14-foot wind tunnel, and the flagged symbol 


On— — r 2 (l .089654 4 In r) 
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represents llto latter result corrected for wall 
interference. Tt can be seen that the latter result 
is very close to the theoretical curves computed 
by means of either equation (81) or (82), but is 
noticeably less than that indicated by equation 
(83). 

Comparison with results for subsonic and super- 
sonic flow. -The remainder of the discussion of 
the aerodynamic properties of cone-cylinders is 
concerned with an examination of the relation 
between the results given in the preceding section 
for Mach number 1 and vicinity and those avail- 
able from other sources for subsonic and supersonic 
flows. The discussion is of particular interest 
because possession of knowledge of the residts for 
Mach numbers near unity greatly facilitates an 
understanding of the process by which the pressure 
distribution changes from the form shown in figure 
15 for incompressible flow to the familiar constant 
pressure associated with supersonic flow. 

In order to be more specific, consider the special 
case of a cone-cylinder having a semiapex angle 
of 1/10 radian. Figure 21 shows the pressure 
distribution on the conical part of such a body for 
the} highest Mach number for purely subsonic 
flow (namely 0), for Mach number 1 , and 



FirjT RK 21. -Pressure distributions on a cone-cvlinder at 
the upper and lower critical Mach numbers and at Mach 
number 1. 


for the lowest Mach number for purely supersonic 
flow. The hitter Mach number is very nearly 
1.055. The results for zero Mach number are 
those given by Fracnkel in reference 31 and shown 
previously herein in figure 15. The results for 
Mach number 1 are those indicated by equation 
(72) for the forward third of the cone and by 
equation (65) for the remainder of the cone. The 
results for the upper critical Mach number are 
determined from the numerical results given by 
Oswatitsch and Sjodin in reference 15 and by 
Slum in reference 16 reinterpreted so as to conform 
with the formulation of the equations of transonic 
flow theory given in equations ( 1 ) through (5). 
Although the differences in the numerical results 
are small, such a change is necessary because the 
basic equations used by Oswatitsch and Sjodin 
and by SI ion are not only slightly different from 
those used herein, hut are slightly different from 
each other. The difference is concerned entirely 
with the quantity A/^fy+l) that appears in 
equations (4) and (5). This quantity is replaced 
by 7+1 in the investigation of Shon and by 
(A/a, 2 - ■1)Z T '»I(T +, — «■*), where a* is the critical 
speed of sound, in t lie investigation of Oswatitsch 
and Sjodin. 

A plot of the variation of pressure coefficient C p 
with free-stream Mach number M*, on the surface 
of an infinite cone having a semiapex angle of 1 /10 
radian, as indicated by the theoretical analyses of 
references 15 and 16 reinterpreted in the manner 
described in the preceding paragraph, is shown in 
figure 22 . A second curve labeled C v is included 
on this plot to illustrate the variation with free- 
stream Much number of the critical pressure 
coefficient associated with the occurrence of sonic 
velo ci t y (i . e . , u j Z7«, - ( 1 — 2 ) /[ A 2 (7 + 1 ) ]) on t h c 

surface of the cone. This curve, which is com- 
puted by use of the expression 


M 2 —l 

r — 2 — - 1 

7+1) 


(84) 


does not pass through 0 at free-stream Mach 
number 1 because of the particular manner in 
which the sonic velocity is expressed in the approx- 
imate equations of transonic flow theory. The 
point of intersection of the two curves determines 
the upper critical Mach number and the associated 
value for the surface pressure coefficient that are 
shown in figure 21 . 


SLENDER-BODY THEORY BASED ON APPROXIMATE SOLUTION OF TRANSONIC FLOW EQUATION 27 



F n jure 22. -Variation, of pressure coefficient with Mach 
number on the surface of an infinite cone. 


If the free-stream Mach number is greater 
than the tipper critical, the flow is supersonic 
everywhere and the results summarized in the 
preceding paragraph apply equally to a cone of 
infinite or finite length. If the free-stream 
Much number is less limn the upper critical, hut 
greater than unity, the flow adjacent to the sur- 
face of the cone is subsonic and the flow field as a 
whole is transonic. In this case, the result sum- 
marized in the precluding paragraph may he ex- 
pected to apply only, to a cone of infinite length, 
or to points near the apex of a cone of finite 
length, as in the cafe of a cone-cylinder. It is 
known, in part-icularj that the velocity must he 
sonic at the sohuh(er of a cone-cylinder if the 
flow is subsonic along the surface of the cone. 
Since intermediate pressures would bo anticipated 
at intermediate points along the surface of the 
conical part of a cone-cylinder, tin* two curves 
shown in figure 22 represent the variation with 
Mach number of the upper and lower bounds 
for the values for the pressure coefficient along 
the conical part of a con e-cylinder having a sonii- 
apex angle of 1/10 radian. The smallest Mach 
number for which a value for the surface pressure 
is indicated in figure 22 represents the limit of 
conical flow past an infinite cone and is the lowest 
Mach number for which the how wave is attached 
to cither an infinite or 1 finite cone having a semiapox 
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Fuicre 23. — Variation of pressure coefficient with Mach 
number at various stations along the length of a cone- 
cylinder of thickness ratio 1/10, as indicated by theory. 


angle of 1/10 radian. At still smaller Mach 
numbers, the how wave is detached, or nonex- 
istent if the free-stream Mach number is less than 
unity, and a stagnation point occurs at the apex. 
Tito approximations of slender-body theory are 
not sufficient to provide a solution that exhibits 
the proper behavior at the apex, however, and 
a logarithmic infinity appears at this point instead 
of the proper value associated with a stagnation 
point. 

Further details of the transition from subsonic 
to supersonic flow arc illustrated by the curves 
shown in figure 23 in which the results shown in 
figure 22 are repeated together with additional 
curves for smaller Mach numbers. Several curves 
are shown, each of which represents the variation 
of the pressure coefficient with Mach number at 
a certain station along the length of a cone- 
cylinder having a semiapex angle of 1/10 radian. 
The additional results shown in figure 23 arc* from 
three different sources. The small circles at zero 
Mach number represent the values indicated 
by the solution given in reference 31 by Fraenkcl 
for incompressible flow. The short linos at Mach 
number 1 represent the values indicated by equa- 
tions (65) and (72) of the present theory. The 
dashed lines at subsonic Mach numbers represent 
the variation of pressure coefficient with Mach 
number calculated by use of the following ex- 
pression derived by combination of tin* results of 
Fraenkcl C p .^ for incompressible flow and the 
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influence of Mach number indicated by equation 
(18) with A/? replaced by 1 —MJ: 

C p = — 2r 2 In y I =Mj+ C Piy (85) 

It should be noted that results calculated in this 
way are not of the same quality as the* remainder 
of the results shown in figure 23, since they repre- 
sent a solution of the linearized equations of 
subsonic flow rather than the nonlinear equations 
of transonic flow theory. One important conse- 
quence of this fact that is not apparent from this 
plot is that the value indicated for the pressure 
coefficient at the shoulder is infinite rather than 
the finite value associated with the occurrence 
of sonic velocity. The subsonic results indicated 
on this plot for points removed from the shoulder 
are probably in at least qualitative agreement 
with those that would be indicated by the solution 
of the nonlinear equations of transonic flow theory, 
provided the free-stream Mach number docs not 
approach too closely to unity. It may be noted, 
in particular, that the effect of Mach number 
indicated by equation (85) is simply to shift 
the entire curve representing the pressure distribu- 
tion. A similar result would also occur if the. 
present method of approximation were applied to 
equation (IS) to determine an approximate solu- 
tion of the equations of transonic flow theory. 
The reason is that the entire term containing A 
is lost in the initial differentiation and is, in 
effect, replaced by an unknown constant upon 
integration of the differential equation. It would 
be necessary, however, to include some additional 
refinement in the details of the solution at the 
shoulder before 4 a result would be obtained that 
exhibits the proper behavior in the vicinity of 
this point. 

Additional insight into the nature of the 
variation with Mach number of the pressure 
distribution on a cone-cylinder can be gained by 
examination of the experimental data of Page for 
a cone-cylinder having a semiapex angle of 7°. 
Figure 24 has been included, accordingly, to sum- 
marize the theoretical and experimental results 
for three representative stations along the length 
of the cone. The plain data points represent the 
values measured in a 2-foot wind tunnel. The 
flagged data points represent the values measured 
on the same model in a 14-foot wind tunnel. It 
can he seen that the tw T o sets of results are in good 


agreement at Mach numbers v T ell removed from 
unity, but that substantial differences occur at 
Mach numbers near unity. As noted previously 
in the discussion of these same data, Page attrib- 
utes these differences to wind-tunnel wall inter- 
ference, and proceeds to determine a theoretical 
estimate of the necessary corrections that must be 
applied in order to simulate free-air conditions 
with free-stream Mach number 1 . The solid data 
points represent the results from the 14-foot wind 
tunnel corrected in this way for wall interference. 
Also included in this sketch are a number of solid 
and dashed curves representing theoretical results 
calculated in the same manner as described in 
connection with figure 23. It can bo seen that the 
experimental data support the remarks of the 
preceding paragraph concerning the general nature 
of the transition from subsonic to supersonic flow 
indicated by t lie various theoretical results. These 
data also indicate some of the difficulties en- 
countered in the det ermination of experiment id 
data for free-stream Mach number 1 that are 
essentially free of wind-tunnel w r all interference. 

Although substantial gaps remain for which the 
drag of a cone-cylinder cannot be computed at the 
present time, it is possible to sketch the variation 
of drag with Mach number from the results and 
discussion of the preceding paragraphs. The 
basic framework is provided by three classes of 
results) namely, those for zero Mach number, those 
for Mach numbers near unity, and those for Mach 
numbers greater than the upper critical. At zero 
Mach number the drag is zero in accordance with 
D’Alembert's paradox. At Mach numbers near 
unity the drag is given by equation (81). At 
Mach numbers greater than the upper critical, the 
drag coefficient is the same as for an infinite cone, 
and is equal numerically to the value for C v shown 
in figure 22. These results are illustrated in figure 
25 by solid lines and a small circle. The dashed 
line at Mach numbers less than the upper critical, 
but greater than the Mach number for shock 
attachment, represents the values for the drag 
coefficient of an infinite cone. It follows from 
the discussion of the pressure distribution given in 
the preceding pages that the drag coefficient for a 
cone-cylinder is less than that for an infinite cone. 
The dashed line at subsonic Mach numbers repre- 
sents the values for the drag coefficient that arc 
associated w ith the pressure distribution indicated 
by equation (85). These values are given by the 
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Fkitre 24. Variation of pressure coefficient with Mach number at various stations along tin* length of a conc-eylincIer 
of thickness ratio 0.1225, as indicated by theory and by experiment. 



Flu the 25. -Variation of pressure-drag coefficient with 
Mach number for a cone-cylinder of thickness ratio 1/10, 
as indicated by theory. 


following relation obtained when equations (13) 
and (85) are combined 

r 7J =- 2 T* In 4=37? (86) 

.-VI though it is apparent from the discussion in the 
preceding paragraphs that the latter results are 
not as reliable as the remainder of the results 


shown in figure 25, they are included because they 
display the proper qualitative trend until the Mach 
number approaches very nea^ to unity. There is 
considerable likelihood that the quantitative 
values are too small at Mach numbers much less 
than unity, however, since the pressure distribu- 
tions on which the calculations are based indicate a 
value for C p of negative infinity at the shoulder 
rather than the proper* finite value associated with 
the occurrence of sonic velocity. 

APPLICATION TO THE FRONT HALF OF A PARABOLIC- ARC 
BODY 

As a first application of the present theory for 
flow with free-stream Mach number near 1 around 
a smooth body, consider the case of the front half 
of a parabolic -arc body of revolution of length /and 
diameter d, The expressions for the ordinates /?, 
the cross-section area S \ and the location of the 
points at which S"(x) vanishes are given by equa- 
tions (25), (26), and (27). Attention is confined 
in this section, however, to the part of the body 
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situated between *=0 and j = //2. The pressure 
distribution for this part of the body is independ- 
ent within wide limits, of the shape of the rear 
and is appropriate, in particular, for either the 
ease in which the rear part of the body continues 
as a parabolic-are body or in which it continues as 
a cylinder extending downstream from x=lf 2. 
For the special ease described above, the ordinary 
differential equation for u on the surface of an 
arbitrary body of revolution given by equation 
(53b) reduces to the following form: 


,l r 

~U MJ- 1 -| 


r 47r 


J'„MJ(y+ ljj 

~ 47 tS 1 1 1 

^S" L 


MJ - 1 S'* MJ(y+\)Se c 

7+1) 4tt 4ttj 


(«T) 

As described in connection with the discussion of 
the properties of equation (53) for smooth bodies, 
the calculation of t lie solution of equation (87) is 
started by use of the assumption that u is analytic 
at the singular point where S"(r) vanishes. The 
solution can thus be expanded in a Taylor's series 
from which the values for u can be readily calcu- 
lated for points slightly removed from the singular 
point. Once these values are known, the solution 
for the remainder of the body can be calculated 
numerically by Milne's method. The correspond- 
ing values for the pressure coefficient can be 
calculated directly therefrom with the aid of 
equation (3). 

Comparison with experimental and other theo- 
retical results. — Experimental results for the 
pressure distribution at frcc-stroam Much number 
1 on the front half of parabolic-arc body of 
fineness ratio 6 followed by a cylindrical section 
have been given in reference 19 by Oswatitscli and 
Kcunc. These results are presented graphically 
in figure 20 together with the theoretical pressure 
distribution calculated in the manner described 
above. It can be seen that the agreement be- 
tween the theoretical and experimental results is 
satisfactory for most purposes. 

Attention is called, however, to the fact men- 
tioned previously that the simpler approximation 
of references 19 and 20 in which \ P is replaced by 
t lie constant indicated by equation (50) also leads 
to theoretical results that agree well with the 



Fir;rnK 2G. Pressure distribution on the front half of a 
parabolic-arc body of fineness ratio 0 at Mach number !, 
as indicated by present theory and by experiment. 


experimental pressure distribution shown in figure 
26. The results of such a calculation are shown 
in figure 26 together with the theoretical and 
experimental results from figure 26. Also in- 
cluded is another curve calculated with constant 
\ P evaluated using equation (50) with a z =l 
rather than %. These two values for a 2 are those 
considered originally in references 19 and 20. It 
is shown in these references that both have merit , 
but the value of a 2 =0.50 is recommended because 
it is superior in the treatment of certain important 
details of the flow field. It is also stated that 
according to the similarity rule, the parameter 
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Figure 27. — Pressure distribution on the front half of a 
parabolic-arc body of fineness ratio 6 at Mach number 1, 
as indicated by present theory, by experiment, and by 
other theoretical results. 

a 2 depends on the shape, but not the thickness 
ratio of the body. 

Keune lias modified these conclusions in more 
recent publications (refs. 21, 22, and 23), and has 
suggested, on the basis of considerations relating 
to the mass flow, the selection of such a value for 
a 2 that the value for \ P agrees with half the value 
for (t+1 )< Pzz / U m at the sonic point. It is further 
shown that this procedure leads to values for a 2 
of 0.455 and 0.28 for parabolic-arc bodies having 
thickness ratios r of 0.146 and 0.073. 

Maeder and Thommen have independently 
proposed a similar procedure, but without the 
factor / 2 , and have presented a curve in reference 


26 from which the resulting values for a 2 can bo 
readily determined for any parabolic-arc body 
having a thickness ratio less than 0.25. This 
curve leads to the following values for a 2 , 0.274, 
0.248, and 0.135, for parabolic-arc bodies having 
thickness ratios of J£, 0.146, and 0.073. Although 
these values should be exactly half those indicated 
by Keune for a body of the same thickness ratio, 
it may be observed that the values are actually 
related in a slightly different manner. The pres- 
ent authors have not sought the source of this 
discrepancy since the differences are much smaller 
than those associated with the factor ){ that exists 
between the intended values. The pressure dis- 
tribution on the front half of a parabolic-arc body 
of fineness ratio 6 has been computed using the 
value for a 2 given by Maeder and Thommen, and 
is presented in figure 27. It is evident from 
examination of the various results shown in this 
figure that the procedure adopted for the selection 
of a suitable value for a 2 has a significant effect on 
the resulting pressure distribution. Although 
certain of these procedures lead to results of good 
quality for the pressure distributions on the front 
half of a parabolic-arc body, it should he recalled 
that no selection of a constant for \ P can lead to a 
reasonable pressure distribution on a cone-cylinder. 

Theoretical pressure distributions for the front 
halves of parabolic-arc bodies of revolution of 
fineness ratio 6 y2, 10, 12, and 14 in flows with 
frec-stream Mach number 1 have also been calcu- 
lated by means of the present theory. These 
results, together with those given above for the 
body of fineness ratio 6, are presented in graphical 
form in figure 28. The numerical values for all 
five eases are given in table IV. It should be 
observed Unit only the values for xjl between 0 
and 0.50 arc to be considered in the present dis- 
cussion. This particular family of bodies has been 
selected because experimental pressure distribu- 
tions are available from (lie investigations of 
references 14 and 11. In both cases, tests were 
conducted on bodies that more nearly approximate 
complete bodies than half bodies, but the results 
should be suitable for comparison with the present 
theoretical results because the pressures on the 
front part of such a body do not depend, at frec- 
stream Mach number 1, on the shape of the rear. 
The results of such a comparison are shown in 
figure 29. Although it is evident from the scatter 
of the experimental data that the relative accuracy 
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Ficjurk 28. — Pressure distributions on the front halves of 
a series of parabolic-arc bodies of various fineness ratios 
at Mach number 1, as indicated by present theory. 

of the experimental results diminishes with increas- 
ing fineness ratio, these comparisons show that 
the theoretical and experimental results are in 
essential agreement. 

The results presented in the preceding para- 
graph indicate that the magnitudes of the values 
for C p are strongly influenced by the fineness ratio. 
The similarity rule given by equation (12) in- 
dicates the existence of such an effect, and figure 
30 is included to illustrate quantitatively the 
degree to which the present results conform to the 
similarity rule. In this figure the values for C v 
for each of the bodies are converted to the corre- 
sponding values for a body of fineness ratio 10 by 
use of equation (12). If the present results con- 


formed perfectly to the similarity rule, the results 
for all five bodies would define a single curve. It 
can be seen that the present results when so pre- 
sented do not quite define a single curve, but the 
differences between the various curves arc small. 
This indicates that the present results for the 
pressure distributions on the front halves of a 
series of parabolic-arc bodies of different fineness 
ratios are related, at free-stream Mach number 1, 
in a manner that very nearly conforms with the 
similarity rule. This conclusion is similar to that 
for the cone-cylinder at free-stream Mach number 
] } hut quite different from that for the parabolic- 
arc body at Mach numbers in the vicinity of the 
critical. 

Pressure drag. Once the pressure distribution 
on the front half of a slender parabolic-arc body 
has been calculated, the corresponding value for 
the pressure drag can be found immediately by 
numerical integration of the general relations 
given by equation (13). The necessary calcula- 
tions have been carried out for each of the five 
cases considered in this section, and the resulting 
values for C^=D/[(paal2)U m 2 S m ] at Mach number 
1 are 


l/d 

6 

<>V 2 

10 

12 

14 

C D 

0. 05988 

0. 02903 

0. 020G3 

0. 01-116 

0. 01031 


These values are plotted as a function of r in figure 
31 to illustrate the effect of thickness ratio on C D . 
It can be seen from examination of these values, 
or by inspection of figure 32 in which C D jr 2 is 
plotted as a function of r, that C D is very nearly 
proportional to the square of r as indicated by 
the similarity rule given in equation (14). 

Experimental values for the pressure drag can 
he obtained by a similar integration of equation 
(13) using the measured values for C p given in 
references 11, 14, and 19, and repeated hero in 
figures 29 and 26. Oswatitsch and Keune have 
already given the results of such a calculation in 
reference 19 for the case of the front half of a 
parabolic-arc body of fineness ratio G followed by 
a cylindrical body. Their result that C D equals 
0.058 compares very favorably with the corre- 
sponding value of 0.05988 indicated by the 
present theory. Experimental values for C D have 
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Figure 20. — Pressure distributions? on the front halves of parabolic-arc bodies of various fineness ratios at Mach number 

1, as indicated by present theory and by experiment. 


also been determined for each of the other five 
bodies for which measured pressure distributions 
are presented in figure 29. The results arc in- 
cluded in figures 31 and 32. Attention is directed, 
however, to the fact that the scatter in the meas- 
ured values for C„ for the bodies of fineness ratio 
10, 12, and particularly 14 is such that considerable 
uncertainty exists in the values for C D . An indi- 
cation of the estimated magnitude of the un- 


certainty due to this cause is provided by the 
vertical lines attached to the symbols. Tt can 
be seen that the uncertainties associated with 
scatter in the pressure-distribution data are not 
sufficiently great to account for the result, obvi- 
ously incorrect, that the drag of the body of fine- 
ness ratio 14 is negative. It appears quite possible, 
on the basis of the results for the cone-cylinder 
presented in reference 32 and shown here in 
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Figure 31. “-Variation of pressure-drag coefficient with 
thickness ratio for the front half of a parabolic-arc body 
at Mach number 1, as indicated by present theory, by 
other theoretical results, and by experiment. 


Figure 30. — Pressure distribution on the front half of a 
parabolic-arc body of fineness ratio 10, as obtained by 
direct calculation using the present theory and by use 
of the similarity rule. 

figures 17 and 24, that this property of the experi- 
mental data may be associated with the effects 
of wind-tunnel wall interference. The results in 
figures 31 and 32 are in agreement with the fact 
that the effects of wind-tunnel wall interference 
at frcc-stream Mach number 1 are of greater 
importance, relative to the quantities measured, 
for thin bodies than for thick bodies. 

It is of interest at this point to compare the 
experimental and theoretical values for C D given 
in the preceding paragraph with those indicated 
by the linearized theory of sonic flow given in 
references 19 and 20. This theory leads to a simple 
rule for the calculation of the drag of half bodies 
that states that the drag docs not depend on the 
value selected for a 2 in equation (50) and is equal 
to exactly half of the value indicated by linearized 
supersonic theory. This rule, when applied to 
the front half of a parabolic-arc body of thickness 
ratio r, results in the following simple expression 
for C D : 

f?z>— (7/3) r 2 (88) 


the values indicated by equation (88) are some- 
what larger than cither the experimental values 
or the calculated values indicated by the present 
theory. Before leaving this topic, it should be 
mentioned that Keunc states briefly (ref. 23) that 
a more elaborate analysis of this case in which 
the quantity a 2 is replaced by a function of x , 
rather than by a constant, leads to a value for 
C D that is about 10 percent smaller than indicated 
by equation (88). This result corresponds very 
closely with the results of the present calculations. 

The approximate theory of Miles (ref. 39) 
described briefly here in the section on cone- 


co 



Curves illustrating the variation of C D with r 
indicated by this result have been included in 
figures 31 and 32, from which it can he seen that 


Figure 32. — Variation of Cj>/t 2 with thickness ratio for 
the front half of a parabolic-arc body at Mach number 
1, as indicated by present theory, by other theoretical 
results, and by experiment. 
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cylinders also contains a general expression for 
the calculation of the pressure drag of the forepart 
of a smooth slender body of revolution in a flow 
with free-stream Mach number 1. Miles’ expres- 
sion, which appears to lead, in general, to very 
small values for the drag of smooth bodies, reduces 
to the following expression when applied to the 
front half of a parabolic-arc body: 

C Y i > = 0.0735r 2 (89) 

This result is also plotted in figures 31 and 32 
from which it is immediately apparent that these 
values for C D are much smaller than either the 
experimental values or the calculated values 
indicated by the present theory or by the theory 
of Oswatitsch and Keune. 

APPLICATION TO COMPLETE PARABOLIC-ARC BODY 

The calculation of the pressure distribution 
along the entire length of a complete parabolic-arc 
body for free-stream Mach number 1 represents an 
interesting application of the procedures described 
in the preceding pages. The pressure distribution 
on the front half of such a body is, of course, the 
same as discussed at length in the preceding sec- 
tion. Although the procedures can be used to 
calculate the pressures at points somewhat rear- 
ward of the middle of the body, they cannot be 
used to calculate the pressures along the entire 
length of the body. The reason is that the velocity 
reaches a maximum at a point somewhat rearward 
of the middle, and then decreases continuously 
along the remainder of the body, thereby violating 
the condition imposed by equation (46) that du/bx 
or X P be positive. The flow, moreover, is subsonic 
in the vicinity of the rear tip. 

The pressure distribution along the entire length 
of a parabolic-arc body can, nevertheless, be calcu- 
lated by considering the solution in sections and 
joining together the various results in the same 
manner as described for t lie cone-cylinder. This 
procedure is quite analogous to the procedures 
commonly employed in many other theoretical in- 
vestigations of transonic flows in which the condi- 
tions in the subsonic region and a restricted part of 
the adjacent supersonic region arc computed using 
hodograph or relaxation methods and the solution 
for the remainder of the supersonic region is com- 
puted by the method of characteristics. The fol- 
lowing discussion is organized accordingly into 


three parts concerned, respectively, with transition 
through sonic velocity on the forebody, supersonic 
flow on the afterbod}", and subsonic flow in the 
vicinity of the rear tip. 

Transition through sonic velocity on the fore- 
body. — A principal feature of the flow on the 
forward part of a body is the smooth transition 
from subsonic velocities at the nose to supersonic 
velocities near the middle of the body. It is amply 
demonstrated in the preceding section that equa- 
tion (87) can be used together with equation (3) to 
calculate the pressure distribution on the front half 
of a parabolic-are body of revolution. The same 
procedures can also be applied to calculate the 
pressures at points somewhat rearward of the 
middle. Such calculations have been made for 
each of the five parabolic-arc bodies of fineness 
ratio 6, G y2. 10, 12, and 14 for which experimental 
results are available from references 11 and 14. 
The calculated results are summarized in table IV. 
The pressure distribution cannot be calculated in 
this way for points along the entire rear half of t he 
body since, as can be observed by examination of 
equation (87), a discontinuity in dufdx, and hence 
the slope of the pressure distribution curve, occurs 
at the point on the rear of the body where S"(x) 
vanishes. This equation shows, moreover, that 
du/dx is zero just upstream and infinite just down- 
stream of this point. The property is symptomatic 
of the breakdown of the approximation furnished 
by equation (87) and signals that the use of this 
relation should be terminated at some point more 
forward along the length of the body. As a conse- 
quence, the results enumerated in table IV for four 
of the five bodies considered are terminated well 
forward of the point on the rear of the body where 
S"(x) vanishes. The results for the fifth body, 
namely that of fineness ratio 10, are given for all 
points back to that at which S"(x) vanishes in 
order to illustrate certain features in the discussion 
of the supersonic flow on the afterbody. The 
results for the body of fineness ratio 10 arc shown 
graphically in figure 33. 

Supersonic flow on the afterbody. — The calcula- 
tion of the pressure distribution on the portion of 
the remainder of the body along which the local 
velocity is supersonic can be continued by appli- 
cation of equation (42) rather than equation (87). 
Inasmuch as these two equations are based on 
somewhat different approximations, the two sets 
of results cannot be expected to continue analyti- 
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Figure 33. — Pressure distribution on the forward part of 
a parabolic-arc body of revolution at Mach number 1 , 
as indicated by equation (87). 


cally into each other at the point at which the t wo 
sets of results are joined. The determination of 
the point at which both C P and dCJdx match can 
be accomplished in a simple manner once a series 
of values for u for various stations along the length 
of the body lias been calculated by numerical inte- 
gration of equation (87). Substitution of these 
values into equations (42) and (87) leads to two 
sets of values for dujdx , each of which defines a 
single curve when plotted as a function of xjl as 
illustrated in figure 34 for the special case of the 



Figure 34. — Illustration of graphical technique for select- 
ing the point at which the solution obtained by equation 
(42) is to be joined to that obtained by equation (87). 


parabolic-arc body of fineness ratio 10. As can be 
seen, the two curves intersect at a certain point, 
in this case situated at a value for xjl of about 
0.542. Since the two values for u and dujdx and 
hence C p and dCJdx arc the same at the point of 
intersection whether calculated using equations 
(42) or (87), this point determines the desired loca- 


tion of the point at which t lie two solutions are to 
be joined. The procedure lias been followed in the 
present calculations except for the slight simplifi- 
cation introduced by joining the solutions at the 
nearest point for which values for C p are given in 
table IV. This simplification results in a slight 
mismatching of the values for dC p jdx , but the 
effects arc small since it follows from the fact that 
these values for C v are given for every 2 percent of 
the body length that the point at which the two 
solutions are joined is at most only 1 percent of 
the body length from the point at which the two 
values for dC p jdx are exactly the same. 

Once the point at which the solutions are to be 
joined and the associated value for u at this point 
are determined, the calculation of the pressure dis- 
tribution on the remainder of the region of super- 
sonic flow on the afterbody can proceed in a 
straight forward manner by application of any of 
the many standard methods available for the nu- 
merical integration of equation (42). Table V pre- 
sents the results of such calculations performed by 
the method of Milne (ref. 9) described in a previ- 
ous section of this report, for each of the five 
parabolic-arc bodies of fineness ratio 6, 6 y2, 10, 
12, and 14. The calculations for the body of fine- 
ness ratio 10 were carried out for points situated 
upstream as well as downstream of the point at 
54.2 percent of the body length at which the solu- 
tions were joined in order to investigate the degree 
to which the two approximate solutions overlap. 
The results are shown graphically in figure 35. It 



Figure 35. — Pressure distribution on a parabolic-arc 
body of revolution at Mach number I as indicated by 
equation (87) and by equation (42). 


can be seen that the two sets of results agree very 
well over a substantial fraction of the length of the 
body. They indicate that the selection of the 
point at which the solutions are joined is not crifri- 



SLEXDER-BODY THEORY BASED OX APPROXIMATE SOLTTTIOX OF TRAXSONIC FLOW EQUATION 37 


cal, and that the approximation afforded by use of 
cither equation (42) or (87) is ver}^ nearly the 
same provided the values for neither u nor bujbx 
approach too close to zero. 

It can be seen from the. results presented above 
that the region of supersonic flow docs not extend 
all the way to the rear tip of the body. Instead, 
the flow decelerates and reaches sonic velocity 
at a point somewhat forward of the rear tip, for 
example at 8G.3 percent of the body length. 
The procedures described above terminate at 
this point since u is zero and the quantity 
ln(Af» 3 — 1 +&?/) in equation (42) leads, at Jl/*— 1, 
to a logarithmic infinity in dujdx . 

The termination of the supersonic region for- 
ward of the rear tip of a parabolic-arc body and 
the great increase in the value of dujdx as the 
sonic point is approached arc similar to the be- 
havior indicated both by the present calculations 
for Mach numbers greater than unity and by 
results calculated by the method of characteristics 
in reference 42. Most of the examples presented 
arc for Mach numbers considerably greater than 
unity, but one example is included in this reference 
in which an approximate solution is determined 
for a parabolic-arc body of fineness ratio G in a 
flow with free-stream Mach number 1 . Inasmuch 
as the method of characteristics is appropriate 
only in a region where the flow is supersonic, it is 
necessary in this particular example to introduce 
some other procedure for the determination of 
the solution for the region of subsonic flow that 
exists in the vicinity of the forward part of the 
body. It is important that a good approximation 
be used since the results not only determine the 
values for the subsonic region, but also affect the 
values calculated by the supersonic region as well. 
The reason is that the initial values required to 
start the calculations with the method of charac- 
teristics at the upstream boundary of the super- 
sonic region must he provided by the solution for 
the subsonic region. The particular method 
used by Oswatitsch was the linearized theory for 
sonic flow described in references 19 and 20 which 
follows from the equations giver herein upon 
replacing A P in equation (47) by the constant 
value indicated by equation (50). The results 
given by Oswatitsch arc shown in figure 3G 
together with the corresponding results indicated 
by the present calculations for tin* same body. 
It can be seen that the two theoretical results arc 



Figure 3G.— Pressure distribution on a parabolic-arc 
body of revolution at Mach number 1 as indicated by 
present theory and by other theoretical results. 


in essential agreement with regard to the details 
of the behavior of the solution in the vicinity of 
the rear sonic point. The slight discrepancies 
between the two results over the forward mid 
middle portions of the body arc associated with 
the selection of an appropriate value for \ P , and 
would be less if the calculations were repeated 
using a smaller value for \ P as advocated by 
Keune in references 21 , 22, and 23 and by Maeder 
and Thornmcn in reference 2G, as discussed in 
connection with figures 2G and 27. 

Subsonic flow near rear tip. -The results pre- 
sented in the preceding section clearly indicate 
that the flow along the surface of a parabolic-arc 
body of revolution decelerates smoothly and 
continuously from a maximum supersonic velocity 
at about two-thirds of the body length to sonic 
velocity at a point definitely forward of the rear 
tip. It is known, furthermore, that the flow in 
the immediate vicinity of the rear tip must be 
subsonic. Little more has been established pre- 
viously regarding the nature of the solution in this 
region. It is evident, however, that two possi- 
bilities exist regarding the nature of the transition 
from supersonic to subsonic flow along the body 
surface. The transition may be accomplished in 
a discontinuous manner involving one or more 
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shock waves that extend to the body surface, as 
is familiar in two-dimensional flows at Mach 
numbers somewhat greater than the critical (see 
ref, 13 for a theoretical analysis of such cases); 
or it may be accomplished in a continuous manner 
with a smooth deceleration through sonic velocity, 
as is indicated by the first- and second -order 
approximations for supersonic flow past a para- 
bolic-arc body illustrated in figures 11 and 12. It 
should be understood that shock waves are present 
in the flow field in both cases, but do not extend 
(o the surface of the body in the ease of smooth 
deceleration through sonic velocity. 

It is assumed in the present investigation that 
the transition from supersonic to subsonic flow 
along the surface of the body is of a continuous 
nature and that the results calculated in the 
preceding section are applicable right up to the 
sonic point. It is assumed, furthermore, that 
the results for the subsonic flow along the remain- 
der of the body can bo calculated by means of 
equation (29). This equation is useful for the 
analysis of subsonic flows, but cannot be applied 
to supersonic flows because of the presence of the 
term ln(l — MJ—ku). This term leads, more- 
over, to a logarithmically infinite value of dufdx 
at the sonic point. It follows from these prop- 
erties that the two sets of values for the ordinate 
u and the slope dufdx calculated by use of equa- 
tions (29) and (42) match when joined at the 
sonic point. The two sets of results do not 
overlap, as in the cases described in the preceding 
section, since the regions of applicability of 
equations (29) and (42) both terminate abruptly 
at the sonic point. It. should be remarked 
before proceeding to the discussion of the detailed 
results of the flow at the rear of the body that 
the logarithmically infinite deceleration at the 
sonic point is not altogether unreasonable phys- 
ically since it may be that such a behavior 
represents the vanishing influence of a shock 
wave that exists in the flow field but does not 
quite extend to the body surface. The experi- 
mental investigation of such a detail would be a 
difficult task because the presence of the boundary 
layer would probably modify the phenomenon 
to such an extent that considerable uncertainty 
would be introduced into the interpretation of 
the results. 

The pressure distribution on the rear of each of 
the five parabolic-arc bodies considered in the pre- 


ceding sections has been calculated by use of equa- 
tion (29) in the manner described above. The re- 
sults were determined again by application of 
Milne's method. The only part of the calculation 
that remains to be described is the detail by which 
the starting value was determined so that the solu- 
tion of equation (29) would indicate the same 
sonic point as the solution of equation (49). It is 
immediately evident that a special treatment of 
some sort is^neccssary because the occurrence of 
infinite values for dufdx at the sonic point imposes 
a definite limit on the radius of convergence of a 
power series expansion of the solution such as em- 
ployed herein in the previous applications of 
Milne's method. 

The first step in the procedures that were used 
in the present calculations is to increase the pre- 
cision with which the sonic point is located. This 
is desirable because the results for the supersonic 
region summarized in table Y were calculated for 
intervals of 2 percent of the body length, and 
thereby leave a possible uncertainty of nearly 
equal magnitude in the precise location of the 
sonic point. The desired increase in precision 
could, of course, be attainable by introducing 
smaller intervals in the application of Milne’s 
method, but this is rather inconvenient because it 
is necessary to repeat the starting procedure or to 
use interpolation formulas every time a smaller 
interval is introduced. Although such a reduc- 
tion in interval size can be more readily accom- 
plished with certain alternative numerical meth- 
ods such as the Runge-Kutta method, the desired 
increase in detail near the sonic point was achieved 
in the present calculations by a careful application 
of the numerical and graphical method of isoclines. 
This procedure has special merit in the present 
application since it is not only capable of locating 
the sonic point when proceeding in the downstream 
direction from the last point in the supersonic 
region for which the results have been calculated 
by Milne’s method, but is also very convenient for 
the continuation of the calculation into the sub- 
sonic region for a sufficient distance downstream 
from the sonic point that the remainder of the 
calculation can again be accomplished with the 
numerical method of Milne. 

The nature of the singularity at the sonic point, 
and the details of the manner in which the solu- 
tions of equations (29) and (42) match when 
joined at this point are illustrated in the series of 
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successive enlargements of the isocline plots shown 
in figure 37 for the body of fineness ratio 10. It 
can be seen that there is an extremely small region 
in which the values for du/dx arc large and that 
the values for du/dx are remarkably similar at 
short distances upstream and downstream from 
the sonic point. The method of isoclines was used 
to calculate the solution for the portion of the 
body that extends from about 85 to 89 percent of 
the body length. The results for the remainder 
of the body were calculated using the method of 
Milne with intervals of 1 percent of the body 
length. The numerical results are summarized in 
table VI, and a plot of the pressure distribution 
for the complete parabolic-arc body of fineness 
ratio 10 is shown in figure 38. It can be seen from 
this plot that the local effects associated with the 
infinite deceleration at the sonic point are of such 
small magnitude and extent that they cannot be 
represented on a plot of this scale. The corre- 
sponding results for the rear parts of the other 
parabolic-arc bodies of fineness ratio 6, 6 ^2, 12, 
and 14 have also been computed in the same way. 
The results of these calculations are also summa- 
rized in table VI. 

Summary and discussion of results. — The calcu- 
lations described in the preceding sections provide 
a set of values for the pressure coefficient C p at 
free-stream Mach number 1 for points spaced at 
intervals no greater than 2 percent of the body 
length along the surface of each of five parabolic- 
arc bodies having fineness ratios of 6, 6^/2, 10, 12, 
and 14. These results are shown graphically in 
figure 39 together with the experimental results 
from references 11 and 14 for parabolic-arc bodies 
truncated at about 5/6 of the length, as indicated 
by the dashed lines. It can be seen that the the- 
oretical and experimental results arc in essential 
agreement, and that the principal discrepancies 
arc confined in each case to a small region in the 
immediate vicinity of the base of the body. This 
discrepancy is undoubtedly associated, for the 
bodies of fineness ratio 6 and 6y2, with the occur- 
rence of a shock wave that must be detached from 
the corner because the local Mach number is too 
small. The situation is not so clear for the other 
bodies, although it would appear from inspection 
of the data shown in figure 5 that the phenomena 
involved are not too different from those for the 
bodies of fineness ratio 6 and 6^2. 

Comparison of the calculated results for the 


parabolic-arc bodies of different fineness ratio 
reveals that the values for the pressure coefficient 
C p vary with thickness ratio r in a manner that is 
very nearly in accordance with the similarity rule 
given by equation (12). This fact is illustrated 
in figure 40 in which are plotted the pressure dis- 
tribution for a parabolic-arc body of fineness ratio 
10 as determined by direct calculation, and by use 
of the similarity rule together with the calculations 
for the bodies of fineness ratio 6, 6^/2, 12, and 14. 
The five curves would coalesce to form a single 
curve if the results of the present calculations were 
in perfect agreement with the similarity rule. It 
can be seen that the five curves do not coalesce 
perfectly onto one line, but the deviations there- 
from are small for the range of thickness ratios 
considered. This result is very similar to that 
stated previously in the discussion of the cone- 
cylinder at free-stream Mach number 1, but quite 
different from the behavior of the present results 
for Mach numbers in the vicinity of the critical. 

Pressure drag. — A value for the pressure drag 
for each of the five parabolic-arc bodies for which 
the pressure distribution is given in the preceding 
section has been calculated by numerical integra- 
tion of equation (13). The resulting values for 
C D =D/[(p ao /2)UJS m ] at Mach number 1 are 


lid 

6 

6V2 

10 

12 

14 

Cd 

0. 2485 

0. 1288 

0. 09343 

0. 06534 

0. 04675 


These values are plotted as a function of r in figure 
41. The corresponding values of C d /t 2 are plotted 
as a function of r in figure 42 to illustrate that 
C D is very nearfy proportional to the square of r as 
indicated by the similarity rule given by equa- 
tion (14). 

Two additional sets of values for C D for com- 
plete parabolic-arc bodies are also included in 
figures 41 and 42 for purposes of comparison. 
Both of these should be considered as more ap- 
proximate than the present results since they are 
derived through considerations involving the 
linearized theory for sonic flow given by Oswatitsch 
and Keune in references 19 and 20. Straight- 
forward application of that theory to a body that 
is pointed at the nose and is either pointed or 
cylindrical at the stern leads to the simple con- 





Figure 37. -Illustration of graphical technique for obtaining solutions near the rear sonic point. 
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Figure 38. — Pressure distribution on a parabolic-arc 
body of revolution at Mach number 1 as indicated by 
present theory. 


elusion that (he value for C D at free-stream Mach 
number 1 is exactly half that indicated for super- 
sonic Mach numbers by linearized compressible 
flow theory. This result holds, independently of 
the particular value selected for a 2 in equation (50), 
and leads to the following expression for C D for a 
complete parabolic-arc body in a flow with free- 
stream Mach number 1. 

C D = (16/3)t 2 (90) 

Oswatitscli and Keune mention the general con- 
clusion stated above as a purely formal conse- 



(a) l/d—6 (b) i/<Z=6v2 (c) Z/d-10 (d) Z/d=12 (c) l/d=U 

Frc.rRE 39. — P rcs.su rc* distributions on parabolic-arc bodies of various fineness ratios at Mach number I as indicated by 

present theory and by experiment. 
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Figure 40. — Pressure distribution on a parabolic-arc 
body of fineness ratio 10, as obtained by direct calcula- 
tion using the present theory and by use of the similarity 
rule. 



Figure 41.- — Variation of pressure-drag coefficient with 
thickness ratio for a parabolic-arc body at Mach num- 
ber 1, as indicated by present theory and by other 
theoretical results. 
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Figure 42. — Variation of Coh 2 with thickness ratio for a 
parabolic-arc body at Mach number 1, as indicated by 
present theory and by other theoretical results. 


quence of the theory, but suggest that approval of 
the result be withheld for bodies with a pointed 
stern until a check is made to ascertain whether 
the assumption of approximately constant accel- 
eration would hold near the rear of such bodies. 
It is evident from inspection of the theoretical or 
experimental pressure distributions summarized 
in figure 5 that this condition is not approached, 
even in a qualitative sense, since Zni/dx is positive 
along the forward part of the body, but negative 
along the rear part. The same property is ex- 
hibited by the related results calculated by the 
method of characteristics that are shown in 
figure 36, as well as by similar results given in 
reference 43 for additional examples. It is shown, 
furthermore (refs. 25 and 26), that even the results 
calculated by use of the linearized theory for sonic 
flow display negative values for bu/bx over the 
rear of a parabolic-arc body in spite of the fact 
that dw/dx is replaced by a positive constant in 
the derivation. Maeder and Thommen have 
nevertheless advocated in references 25 and 26, 
that equation (88) be used to calculate the pressure 
drag of a parabolic-arc body of revolution in a 
flow with free-stream Mach number 1. It can be 
seen that the values for Cd indicated by equation 
(90) are much smaller than those indicated by the 
present calculations. 

The second additional set of values for C D 
included in figures 41 and 42 has been calculated 
by a simple rule proposed in reference 43. It 
states that the value for C D for a body of revolu- 
tion in a flow with free-stream Mach number 1 
is equal to the sum of half the supersonic drag of 
the forebody plus the full value of the supersonic 
drag of the afterbody. The values for the super- 
sonic drag are to be calculated in both instances 
by means of the linearized theory of compressible 
flow, and the line of demarcation between the fore- 
body and afterbody is taken to be the station of 
maximum diameter. This rule is not established 
by a deductive process, but is proposed as a 
generalization, approximate in nature, of observa- 
tions of specific results calculated for three different 
bodies. Application of this rule to the specific case 
of a parabolic-arc body of revolution leads to the 
following expression for C D : 

C D — (25/3) t 2 (91) 
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It can be seen upon inspection of figures 41 and 
42 that the values for C D indicated by equation 
(91) are somewhat smaller than those indicated 
by the more elaborate calculations of the present 
investigation. It should be pointed out that this 
comparison does not fully represent an additional 
test of the rule, however, since the parabolic-arc 
body considered here is one of the three cases 
used by Keune and Oswatitsch in the establish- 
ment of the rule. It appears to the present 
authors, in spite of the moderate success of this 
rule in the three cases considered to date, that the 
final decision regarding the general applicability 
of this rule should be withheld until either similar 
comparisons are made for a much larger class of 
bodies or a better understanding of the funda- 
mental basis for the rule is discovered. 

Comparison with results for subsonic and 
supersonic flow, — The remainder of this report is 
concerned with a discussion of the relation be- 
tween the results calculated for a parabolic-arc 
body of revolution for Mach number 1 and those 
calculated for other Mach numbers. The results 
for both pressure distribution and drag will be 
considered, but attention is confined to the case 
of a body of fineness ratio 10. Although all of 
the results required to make similar comparisons 
for bodies of other fineness ratios have not been 
determined, it is presumed that the relations 
among the results for various Mach numbers will 
remain qualitatively the same. Figure 43 shows 
a comparison of the pressure distribution cal- 
culated for Mach number 1 and those calculated 
for the upper and lower critical Mach numbers. 
It can be seen that the results for Mach number 1 
are qualitatively similar to those for the lower 
critical Mach number along the forward part of 
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Figure 43. — Pressure distributions on a parabolic-arc 
body of revolution at the upper and lower critical Mach 
numbers and at Mach number 1. 


the body, and qualitatively similar to those for 
the upper critical Mach number along the rear 
part of the body. The results for various Mach 
numbers are displayed in a second manner in 
figure 44 which shows the variation of the pressure 
coefficient with Mach number for several selected 
points along the surface of the body. Two addi- 
tional curves are included on this plot to show the 
variation with Mach number of the pressure 
coefficients corresponding to the occurrence of 

sonic velocity (i.e., u[U m = ^ ~ at tllC 

midpoint and at the nose of the body. It can 
be seen that the line segments for the surface 
pressures could be connected by reasonably 
smooth curves for points along the front part of 
the body, but not for points along the rear part. 
This situation is similar to that encountered in 
the analysis of transonic flow around airfoils (see, 
e.g., ref. 13) and it is known that the necessary 
adjustment is achieved, in an in viscid flow, in a 
discontinuous manner involving a shock wave 
that extends to the surface of the body. The 
presence of the boundary layer in a viscous flow 
softens the discontinuous change indicated by the 
inviscid theory, but the general behavior is other- 
wise much the same. The actual existence of 
such effects for the special case of a parabolic-arc 
body of revolution is clearly demonstrated experi- 
mentally by the data reported in references 
1 1 and 14. 

The variat ion with Mach number of the pressure 
drag coefficient for a parabolic-arc body of fineness 
ratio 10 is shown in figure 45. The small circles 
represent the values for the drag coefficient of 
the complete body as indicated by the present 
theory. The numerical values are 


Moo 

1. 21 

1. 22 

1. 23 

1. 25 

1. 30 

1. 40 

C D 

0. 1153 

0. 1149 

0. 1146 

0. 1142 

0. 1133 

0. 1122 


The dashed line of figure 45 that is adjacent to 
the small circles represents the values calculated 
by use of the following expression that results 
from both first- and second-order supersonic 
slcnder-body theory. 

C D = (32/3) t 2 (92) 
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Figure 44. — Variation of pressure coefficient with Mach number at various stations along the length of a parabolic-arc 

body of fineness ratio 10, as indicated by present theory. 



Figure 45.- Variation of the pressure-drag coefficient 
with Mach number for a parabolic-arc body of fineness 
ratio 10, as indicated by present theory and by other 
theoretical results. 

It can be scon that the values indicated by the 
present theory are slightly greater than those 
indicated by equation (92). 

A second set of calculated results, indicated 
by small squares, is included in figure 45 to repre- 
sent the contribution to the drag that arises 
from the forces on the front half of the complete 
body. The numerical values are 


Moo 

1. 21 

1. 22 

1. 23 

1. 25 

1. 30 

1. 40 

Co 

0. 04943 

0. 04922 

0. 04904 

0. 04878 

0. 04832 

0. 04783 


The adjacent dashed line in this figure represents 
the values calculated by means of the following 
relation provided by first-order slender-body 
theory 

C D =( 14/3)t 2 (93) 

This result shows once again that only about half 
of the drag rise of a half body occurs at Mach 
numbers less than unity; whereas the previous 
result shows that nearly all of the drag rise of a 
complete parabolic-arc body occurs in the same 
range of Mach numbers. 

A third set of calculated results, indicated by 
small diamonds, is also included in this figure to 
represent the contribution to the drag coefficient 
that arises from the forces on the front three- 
fourths of a parabolic-arc body. The numerical 
results are 
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M ro 

1. 21 

1. 22 

1. 23 

1. 25 

1. 30 

1. 40 

Co 

0. 09508 

0. 09455 

0. 0940G 

i 

0. 09320 

0. 09147 

0. 08898 


The adjacent dashed lino represents the values 
indicated by the first-order slender-body theory 
which have been calculated with the aid of the 
following expression : 



VAf. 2 — 1 


■) 


(94) 


Although it is clear from this expression that the 
drag is not independent of the free-stream Mach 
number, as in the case of the complete body and 
also the half body, it can be seen from figure 45 
that the effect is very small. 

CONCLUDING REMARKS 

It has been amply demonstrated in the preceding 
pages that the approximate method for the solu- 
tion of the nonlinear equations of transonic flow 
theory described in reference 4 can also be ap- 
plied to axisymmetric flows. As in the case of two- 
dimensional flows, the results are shown, by com- 
parison with experimental and other theoretical 
results, to be of almost surprising accuracy con- 
sidering both the small perturbation and slender- 
body approximations inherent in the fundamental 
equations of the theory and the novel nature of 
the procedures used to obtain a simple approxi- 
mate solution of these equations. 

The results of this investigation are of interest 
not only because of the frequent use of a body of 
revolution in aeronautical design, but also because 
of the central role of the body of revolution in 
applications of the transonic area and equivalence 
rules. The best known of these ndes is probably 
that for wave drag discovered experimentally by 
ATI) it comb and announced in reference 1; but 
recent developments (ref. 3 for a resume) enable 
the rcad}^ calculation of not only the drag, hut also 
the pressure distribution and other aerodynamic 
properties of slender wings, bodies, wing-body 
combinations, etc., provided that the pressure 
distribution is known for the equivalent nonlifting 
body of revolution having the same longitudinal 
distribution of cross-section area. Thus the 
utility of these rules is greatly increased by the 


availability of the methods described in (bis paper 
for the calculation of the pressure distribution on 
the equivalent body. Although no examples of 
such an application arc given in the present paper, 
the results for a simple case involving a thin elliptic 
cone that resembles a delta wing are given in 
reference 5. 

While the development of a reliable theory for 
the prediction of the aerodynamic properties of 
bodies of revolution in transonic flow is a subject 
primarily of interest to acrodynamicists, it would 
seem to the authors that the method of approxima- 
tion and some of the results for specific applica- 
tions might also be of interest to mathematicians 
concerned with the solution of nonlinear partial 
differential equations. It should be remarked in 
this connection, that the method of approximation 
is only described briefly and from a simple heuristic 
point of view, in the present paper. Those inter- 
ested in a further discussion of the method from 
a number of alternative points of view are referred 
to the appendix of reference 4. The latter dis- 
cussion is for two-dimensional flow, but the neces- 
sary extension to axisymmetric flow can be ac- 
complished directly without the introduction of 
any new concepts. While it is hoped that the 
discussion of the method of approximation pre- 
sented in reference 4 is illuminating and helps to 
make the procedure more plausible than does the 
brief account presented herein, the underlying 
reason for the success of the method remains only 
incompletely understood from a mathematical 
point of view. On the other hand, the demon- 
strated ability of the method to nearly reproduce 
the results of available, but far more complicated, 
theories and of experiments should be sufficient to 
warrant further mathematical study. 


Ames Research Center 

National Aeronautics and Space Administration 
Moffett Field, Calif., A r ov. 26, 1958 

REFERENCES 

1. Whitcomb, Richard T. : A Study of the Zoro-Lift 

Drag-Rise Charateristics of Wing-Body Combina- 
tions Near the Speed of Sound. NACA Rep. 1273, 
1956. (Supersedes NACA RM L52H08) 

2. Oswatitsch, Klaus, and Keune, Friedrich: Ein Xqui- 

valensatz fur nichtangestellte Flugcl kleiner Span- 
nweitte in schallnaher Stromung. Zeitschrift fur 
Fhigwisscnschaften, 3 Jahrgang, Heft 2, Feb. 1955, 
S. 29-46. 




TECHNICAL REPORT R— 2 — NATIONAL AERONAUTICS AND SPACE ADMINISTRATION 


46 

3. Heaslet, Max. A., and Spreiter, John R.: Three-Di- 

mensional Transonic Flow Theory Applied to 
Slender Wings and Bodies. NACA Rep. 1318, 1957. 
(Supersedes NACA TN 3717) 

4. Spreiter, John R., and Alksne, Alberta Y.: Thin Air- 

foil Theory Based on Approximate Solution of the 
Transonic Flow Equation. NACA TN 3970, 1957. 

5. Spreiter, John R., and Alksne, Alberta Y.: Aerody- 

namics of Wings and Bodies at Mach Number One. 
Proc. Third U.S. National Congress of Applied 
Mechanics, ASME, 1958, pp. 827—835. 

6. Oswatitsch, Klaus: Gas Dynamics (English Version). 

Academic Press, Inc., 1956. 

7. Oswatitsch, Klaus, and Berndt, Sune B.: Aerodynamic 

Similarity at Axisymmetric Transonic Flow Around 
Slender Bodies. KTIt-Aero TN 15, Roy. Inst. 
Tech., Division of Aeronautics, Stockholm, Sweden, 
1950. ' 

8. Ward, G. N.: Linearized Theory of Steady High- 

Speed Flow. Cambridge Univ, Press, 1955. 

9. Milne, William Edmund: Numerical Calculus. Prince- 

ton Univ. Press, 1949. 

10. Van Dyke, Milton D.: Second-Order Slcnder-Body 

Theory — Axisymmetric Flow. NACA TN 4281, 
1958. 

11. Taylor, Robert A., and McDevitt, John B.: Pressure 

Distributions at Transonic Speeds for Parabolic- Arc 
Bodies of Revolution Having Fineness Ratios of 1 0, 
12, and 14. NACA TN 4234, 1958. 

12. Spreiter, John R., and Alksne, Alberta: Theoretical 

Prediction of Pressure Distributions on Nonlifting 
Airfoils at High Subsonic Speeds. NACA Rep. 
1217, 1955. (Supersedes NACA TN 3096) 

13. Spreiter, John R,, Alksne, Alberta Y., and Hyett, B. 

Jeanne; Theoretical Pressure Distributions for 
Several Related Nonlifting Airfoils at High Subsonic 
Speeds. NACA TN 4148, 1958. 

14. Drougge, Georg: Some Measurements on Bodies of 

Revolution at Transonic Speeds. Ninth Interna- 
tional Congress of Applied Mechanics, Univ. of 
Brussels, 1957, pp. 70-77. 

15. Oswatitsch, K. } and Sjodin, L. : Kegelige tJberschall- 

stromung in Schallnahe. Osterr. Ing.-Archiv., vol. 
8, no. 4, 1954, pp. 284-292. 

16. Shen, Yung-chung: Similarity Solution for Transonic 

Flow Past a Cone. OSR TN 56-121, Mar. 1956. 

17. Taylor, G. I., and Maecoll, J. W.: The Air Pressure 

on a Cone Moving at High Speeds. Proc. Roy. Soc. 
of London, seT. A., vol. 139, 1933, p. 278. 

18. Staff of the Computing Section (under the direction of 

Zdenek Kopal) : Tables of Supersonic Flow Around 
Cones. Tech. Rep. no. 1, Center of Analysis, Mass. 
Inst, of Tech., Cambridge, Mass., 1947. 

19. Oswatitsch, K., and Keune, F. : The Flow Around 

Bodies of Revolution at Mach Number 1. Proc. 
Conf. on High-Speed Aeronautics, Polytechnic In- 
stitute of Brooklyn, Brooklyn, N.Y., June 20-22, 
1955, pp. 113-131. 


20. Keune, F.: Bericht tiber eine Naherungstheorie der 

Stromung um Rotationskorper ohne Anstellung bei 
Machzahl Eins. DVL Bericht Nr. 3 V estdeutscher 
Verlag, Koln und Opladen, 1955. 

21. Keune, F.: Uber den Kompressibilitatseinfluss bei und 

nahr Machzahl Eins fur Korper kleiner Streckung 
und schlanke RotationskorpeT. Zs. f. Flugwiss. 4 
Jahr., Heft 1/2, Jan/Feb. 1956, pp. 47-53. 

22. Keune, Friedrich: t!ber cine Naherungstheorie zur 

Berechnung der auftricbsloscn Stromung um Sch- 
lanke Korper bei Schallanstromung. Jahrbuch der 
Wisscnschaft lichen Gesellschaft fur Luftfahrt E.\ . 
(WGL), 1955, pp. 176-186. 

23. Keune, Friedrich: A Contribution to an Approximate 

Calculation of the Flow on Pointed Halfbodics of 
Revolution in the Lower Transonic Range. Ninth 
International Congress of Applied Mechanics, Univ. 
of Brussels, 1957, pp. 25-33. 

24. Maeder, P. F., and Thommcn, H. U.: Some Results 

of Linearized Transonic Flow About Slender Air- 
foils and Bodies of Revolution. Jour. Aero. Sci., 
vol. 23, no. 2, Feb. 1956, pp. 187-188. 

25. Maeder, Paul F.: Solutions to the Linearized Equa- 

tion for Transonic Flows and Their Comparison 
with Experiment. Ninth International Congress of 
Applied Mechanics, Univ. of Brussels, 195/, pp. 
15-24. 

26. Maeder, P. F., and Tliommen, H. U.: Linearized 

Transonic Flow About Slender Bodies of Revolu- 
tion at Zero Incidence. Technical Rep. W T— 25, 
Brown University, July 1957. 

27. Miles, John W.: On Linearized Transonic Flow Theory 

for Slender Bodies. Jour. Aero. Sci., vol. 23, no. 7, 
July 1956, pp. 704-705. 

28. Laitone, E. V.: Experimental Measurement of Incom- 

pressible Flow Along a Cylinder with a Conical 
Nose. Jour. Appl. Pliys., vol. 22, no. 1, Jan. 1951, 
pp. 63-64. 

29. Laitone, E. V.: The Subsonic Flow About a Body of 

Revolution. Quart. Appl. Math., vol. 5, 194/, pp. 
227-231. 

30. Laitone, E. V.: The Linearized Subsonic and Super- 

sonic Flow About Inclined Slender Bodies of Rev- 
olution. Jour. Aero. Sci., vol. 14, 1947, pp. 631-642, 

31. Fraenkel, L. E.: Incompressible Flow Past Quasi- 

Cylindrical Bodies and Some Associated Problems. 
Quart. Jour. Mech. and Appl. Math., vol. XI, pt. 
2, May 1958, pp. 212-222. 

32. Page, William A.: Experimental Study of the Equiv- 

alence of Transonic Flow About Slender Cone- 
Cylinders of Circular and Elliptic Cross Section. 
NACA TN 4233, 1958. 

33. Yoshihara, Hideo: On the Flow’ Over a Cone-Cylinder 

at Mach Number One. WADC Tech. Rep. 52-295, 
Nov. 1952. 

34. Kusukawa, Ken-Ichi: On the Compressible Flow Over 

a Slender Body of Revolution with a Flat Base 
Placed in an Unbounded Fluid and in a Choked 
Wind Tunnel. Jour. Phys. Soc. Japan, vol. 12, no. 
9, Sept. 1957, pp. 1042-1048. 


SLENDER-BODY THEORY BASED ON APPROXIMATE SOLUTION OF TRANSONIC FLOW EQUATION 47 


35. Kusukawa, Ken-Ichi: On the Transonic Flow of a 

Compressible Fluid Past an Axisymmetric Slender 
Body at Zero Incidence. Jour. Phvs. Soc. Japan, 
vol. 12, no. 4, April 1957, pp. 401-410. 

36. Kusukawa, Ken-Ichi: On the Two-Dimensional Com- 

pressible Flow Over a Thin Symmetric Obstacle 
with Sharp Shoulders Placed in an Unbounded 
Fluid and in a Choked Wind Tunnel. Jour. Phys. 
Soc. Japan, vol. 12, no. 9, Sept. 1957, pp. 1031-1041. 

37. Helliwell, J. B., and Mackie, A. G.: Two Dimensional 

Subsonic and Sonic Flow Past Thin Bodies. Jour. 
Fluid Mech., vol. 3, pt. 1, Oct. 1957, pp. 93-109. 

38. Guderley, K, G. : Theorie schallnaher Stromungen. 

Springor-Verlag, 1957. 


39. Miles, John W.: On the Sonic Drag of a Slender Body. 

Jour. Aero. Sci., vol. 23, no. 2, Feb. 1956, pp. 
146-154. 

40. Guderley, G., and Yoshihara, H.: An Axial-Symmetric 

Transonic Flow Pattern. Quart. Appl. Math., vol. 
4, no. 4, April 1951, pp. 333-339. 

41. Guderley, G., and Yoshihara, H.: Axial-Symmetric 

Transonic Flow Patterns. AF Tech. Rep. 5797, 
USAF AMC, Dayton, Ohio, 1949. 

42. Oswatitsch, K.: Die Bcrechnung wirbelfreier achsen- 

symmetrischer Uberschallfelder. Osterr. Ing.- 
Archiv., vol. 10, no. 4, 1956, pp. 359-382. 

43. Keunc, F., and Oswatitsch, K. : On the Influence of 

the Geometry of Slender Bodies of Revolution and 
Delta Wings on Their Drag and Pressure Distribu- 
tion at Transonic Speeds. KTTT-Aero TN 42, 1956. 


TABLE I. — RESULTS OF CALCULATIONS FOR PARABOLIC-ARC BODY OF FINENESS RATIO 10 AT 
SUBCRITICAL MACH NUMBERS (PURELY SUBSONIC FLOWS) 


til 

(a) =0.70 

(b) M co = 0.80 

(c) Mcc =0.85 

(d) A/® =0.90 

u/tfeo 

A/29 

C „ 

U/Uoo 

A/29 

C P 


A/29 


ujUa o 

A/29 

Cp 

0. 0113249 

-0. 125023 

-0.0001566 

0. 213038 

-0. 130509 

-0.0001565 

0. 222809 

-0. 133829 

-0. 0001565 

0. 229450 

-0. 138118 

-0. 0001565 

0- 238027 

. 0313249 

-.093099 

-.0000166 

. 151054 

-.097572 

-.0000165 

. 159998 

-.100641 

-.0000165 

. 166136 

-.104643 

-.0000165 

. 174141 

. 0513249 

-.074023 

-.0000041 

. 115836 

-.078039 

-.0000041 

. 123868 

-.080822 

-.0000041 

. 129434 

— . 0S449O 

-.0000041 

. 136770 

. 0713249 

-.059443 

-.0000014 

. 089485 

-.062978 

-. 0000014 

.096554 

-.065451 

-.0000014 

.101500 

-.068746 

-.0000014 

.108091 

. 0913249 

-.047376 

-.0000006 

.068030 

-.050412 

-.0000006 

.074102 

-.052557 

-.0000006 

. 078392 

-.055449 

-.0000006 

.084175 

. 1113249 

-.037000 

-.0000003 

.049830 

-.039528 

-.0000003 

. 054884 

-.041331 

-.0000003 

.058490 

-.043790 

-.0000003 

. 063409 

, 1313249 

-.027886 

-.0000002 

. 034024 

-.029900 

-.0000002 

. 038052 

-.031350 

-.0000002 

. 040953 

-.033353 

-.0000002 

. 044959 

. 1513249 

-.019773 

-.0000001 

.020095 

-.021274 

-.0000001 

.023096 

-.022365 

-.0000001 

. 025278 

-.023891 

-.0000001 

.028330 

. 1713249 

-.012495 

-.0000001 

.007705 

-. 013486 

-.0000001 

.009688 

-.014213 

-.0000001 

.011142 

-.015244 

-.0000001 

.013204 

.1913249 

-.005932 


-.003380 

-.006422 


-.002401 

-.006785 


-.001676 

-.007305 


-.000634 

.2113249 

0 


-.013333 

0 


-.013333 

0 


-.013333 

0 


-.013333 

. 2313249 

. 005366 


022282 

. 005841 


-. 023231 

.006198 


023947 

. 006726 


— . 025002 

. 2513249 

. 010216 

.0000000 

030327 

.011147 

.0000000 

-. 032188 

.011854 

.0000000 

-.033602 

.012911 

.0000000 

-.035717 

. 2713249 

. 014589 
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037546 

.015954 

. 0000000 

-.040274 

.016999 

.0000000 

-.042365 

. 018584 
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045534 

. 2913249 

.018517 
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-.044001 

.020290 
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-.047548 

.021659 
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-.050286 

. 023763 
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-.054493 

.3113249 

. 022023 
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049743 

.024179 
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-.054053 

.025854 
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-.057404 

. 028462 
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-.062620 

. 3313249 

.025130 

.0000000 

-.054813 

. 027637 
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-. 059825 
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063749 
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- . 069932 
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-.064894 
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.0000000 
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.0000000 
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. 038228 
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-. 087035 

.4113249 
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-. 068963 
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-.091121 
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. 039982 
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-.080344 

.043078 
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-.086536 

. 048228 
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. 036764 
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073659 

. 040703 

.0000000 

081537 

.043871 
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-.087873 
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. 049579 
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. 5313249 

. 036698 
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073554 

. 040629 
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-.081414 

. 043789 

.0000000 

-.087736 

. 049064 
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-.098286 

. 5513249 

.036018 
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-. 072457 

. 039859 
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i 

-.080139 

i 

.042943 

.0000000 

086307 

.048069 

-.0000000 

-.096559 
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TABLE I. RESULTS OF CALCULATIONS FOR PARABOLIC-ARC BODY OF FINENESS RATIO 10 AT 
SUBCRITICAL MACII NUMBERS (PURELY SUBSONIC FLOWS)— Concluded 



(e) A/co =0.92 

(f) A/® =0.93 

(g) A/qo =0.934 

(h) A/oo =0.935 


U/Uoi 

A/29 

Cp 

U/Uco 

A/29 

Cp 

U/Uco 

A/29 

Cp 

UjUco 

A/29 

Cp 

0. 0113249 

-0.140235 

-0. 0001 565 

0. 242261 

-0. 141412 

-0.0001565 

0. 244616 

-0. 141909 

-0. 0001565 

0. 245610 

-0. 142036 

-0. 0001565 

0. 245863 

. 0313249 

-. 106637 

-.0000165 

. 178129 

107751 

-. 0000165 

.180357 

108222 

-.0000165 

. 181300 

-.108344 

-.0000165 

. 181543 

. 0513249 

— . 086335 

-.0000041 

. 140461 

-.087372 

-.0000041 

. 142534 

- 087812 

-.0000041 

. 143414 

-.087924 

-.0000041 

. 143639 

.0713249 

-.070422 

-.0000014 

.111443 

-.071370 

-.0000014 

.113337 

-.071773 

-.0000014 

. 114144 

-.071876 

-.0000014 

. 114351 

.0913249 

-.056936 

-.0000006 

. 087150 

-.057782 

-.0000006 

.088843 

-.058144 

-.0000006 

. 089566 

-. 058237 

-.0000006 

. 089752 

.1113249 

-. 045071 

-.0000003 

. 065972 

-. 04580G 

-.0000003 

.067440 

-.046121 

-.0000003 

.068070 

-. 046202 

-.0000003 

. 068233 

,1313249 

-.034411 

-.0000002 

. 047074 

-. 035022 

-.0000002 

. 048297 

-.035285 

-.0000002 

. 048824 

- 035353 

-.0000002 

. 048959 

. 1513249 

-.024708 

-.0000001 

. 029963 

025184 

-.0000001 

.030916 

-.025390 

-.0000001 

.031329 

-.025444 

-.0000001 

.031436 

. 1713249 

-.015804 

-.0000001 

.014324 

-016134 

-.0000001 

.014984 

-.016278 

-.0000001 

. 015271 

-.016315 

-.0000001 

.015346 

, 1913249 

007593 


-.000059 

-.007764 


. 000284 

-.007839 


.000434 

-. 007859 


.000473 

.2113249 

0 


-.013333 

0 


-.013333 

0 


-.013333 

0 


-.013333 

. 2313249 

. 007027 


-.025605 

. 007212 


-. 025974 

.007294 


-.026138 

. 007316 


— . 026181 

. 2513249 

. 013527 

.0000000 

-.036948 

.013909 

.ooooooo 

-.037713 

. 014082 

.0000002 

-. 038058 

.014127 

.ooooooo 

-.038148 

. 2713249 

. 019525 

.0000000 

-.047416 

.020119 

.ooooooo 

-.048605 

. 020390 

.ooooooo 

-.049147 

. 020461 

.ooooooo 

049290 

. 2913249 

. 025037 

.0000000 

057042 

.025858 

.ooooooo 

-.058683 

. 020237 

.ooooooo 

-.059441 

. 026338 

.ooooooo 

-.059642 

.3113249 

.030075 

.0000000 

-.065847 

.031135 

.ooooooo 

-.067967 

. 031634 

.ooooooo 

068963 

. 031767 

.ooooooo 

-.069229 

. 33T3249 

. 034643 

.0000000 

-.073838 

.035956 

.ooooooo 

-. 076463 

. 036585 

.ooooooo 

-.077721 

. 036754 

. ooooooo 

-.078061 

. 3513249 

. 038739 

.0000000 

-.081014 

.040315 

.ooooooo 

084167 

. 041088 

.ooooooo 

085713 

.041299 

.ooooooo 

-.086135 

.3713249 

. 042358 

.0000000 

-.087364 

,044206 

.ooooooo 

-.091061 

. 045137 

-.ooooooo 

092923 

.045395 

-.ooooooo 

-.093440 

.3913249 

. 045491 

. ooooooo 

-.092871 

. 047613 

-.ooooooo 

097116 

. 048718 

-.ooooooo 

-.099325 

.049031 1 

-.ooooooo 

-.099951 

.4113249 

.048126 

.0000000 

- 097511 

,050517 

-.ooooooo 

-. 102292 

. 051809 

-.ooooooo 

-.104875 

.052186 

-.ooooooo 

— . 105630 

. 4313249 

.050251 

.0000000 

101256 

. 052890 

.ooooooo 

106536 

. 0543S0 

-.ooooooo 

-.109515 

.054833 1 

-.ooooooo 

-.110421 

. 4513249 

.051850 

.0000000 

-. 104080 

. 054703 

.ooooooo 

109786 

. 056390 

-.ooooooo 

-.113159 

. 056933 

-.0000001 

-.114245 

. 4713249 

. 052913 

.0000000 

-.105957 

. 055922 

.ooooooo 

111977 

. 057779 

-.ooooooo 

-.115690 

. 05S424 

-.0000001 

-.116979 

. 4913249 

. 053428 

.0000000 

-.106868 

. 056519 

.0000001 

113051 

. 058478 

.0000001 

-.116968 

. 059206 

-.0000003 

-.118424 

.5113249 

. 053392 

.0000000 

106804 

, 056477 

.0000001 

112974 

. 058428 

.0000005 

116876 

.059153 

.0000010 

-.118326 

. 5313249 

. 052803 

. ooooooo 

-.105764 

. 055796 

.ooooooo 

-.111749 

. 057633 

.0000004 

-.115424 

. 058261 

.0000021 

-.116680 

. 5513249 

.051669 

-.0000000 

- 103760 

. 054548 

-.0000019 

109518 

. 056158 

-.0000005 

-. 112737 

. 056692 

-.0000024 

113805 
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TABLE III.- -RESULTS OF CALCULATIONS FOR PARABOLIC-ARC BODY OF FINENESS RATIO 10 AT 
MACH NUMBERS ABOVE TIIE UPPER CRITICAL (PURELY SUPERSONIC FLOWS) 


xjl 

(a) M<n —1.21 

(b) A/oo =1.22 

(C) Moo -1.23 


u(U*> 

d/29 

Cp 

It/ L OO 

A/29 


ll/V oo 

A/29 

Cp 

0. 01 13249 

-0. 122523 

-0. 0000057 

0. 206838 

-0.119625 

-0. 0000021 

0. 201041 

-0. 117275 

-0.0000010 

0. 196342 

.0313249 

-.110329 

-.0000011 

. 185513 

108363 

-.0000006 

. 181582 

-.106645 

-.0000003 

.178144 

,0513249 

-.099655 

-.0000004 

. 167101 

-.098195 

-.0000002 

. 164181 I 

-.096881 

-.0000002 

. 161552 

.C713249 

-.089884 

-.0000002 

. 150365 

-.088770 

-.0000001 

. 148139 

-. 087752 

-.0000001 

. 146102 

.0913249 

-. 080765 

-.0000001 

, 134808 

-.079914 

-.0000001 

. 133106 

-. 079127 

-.0000001 

. 131532 

. 1113249 

-.072163 

-.0000001 

. 120154 

-.071520 

-.ooooooo 

. 118870 

-.070922 

-.ooooooo 

.117674 

.1313249 

- 063990 

-.0000000 

. 106232 

-.063520 

-.ooooooo 

. 105292 

-.063080 

-.ooooooo 

. 104412 

. 1513249 

-.056188 

-.0000000 

.092923 

-.055863 

-.ooooooo 

. 092273 

-.055557 

-.ooooooo 

.091662 

. 1713249 

048715 

-.0000000 

. 080145 

-.048514 

-.ooooooo 

. 079743 

-. 048324 

- . ooooooo 

.079364 

. 1913249 

-.041.540 


.067835 

-.041447 


. 067648 

-.041358 


. 067471 

. 2113249 

-.034641 


. 055949 

-.034641 


. 055949 

-.034641 


.055949 

.2313249 

028001 


. 044451 

-.028082 


.044615 

-.028160 


.044771 
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TABLE III. -RESULTS OF CALCULATIONS FOR PARABOLIC-ARC BODY OF FINENESS RATIO 10 AT 
MACH NUMBERS ABOVE THE UPPER CRITICAL (PURELY SUPERSONIC FLOWS) — Concluded 
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TABLE V.— RESULTS OF CALCULATIONS FOR PARABOLIC-ARC BODIES OF VARIOUS FINENESS RATIOS AT = 1 - EQUATION (42) 
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